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PREFACE 


This  report  provides  a  description  of  the  current  status  of  a  millimeter  wave 
propagation  model  based  on  geophysical  input  data.  The  model  is  being  developed  by 
the  National  Telecommunications  and  Information  Administration,  Institute  for  Tele¬ 
communication  Sciences,  to  aid  in  the  development  of  millimeter  wave  telecommunica¬ 
tion  syscems. 

The  laboratory  measurements  that  have  been  used  to  verify,  and  at  times,  drive 
the  model,  have  been  supported  by  the  Army  Research  Office,  under  Contracts 
ARO  i 01 -83,  6-82,  51-81,  and  42-80. 
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AN  ATMOSPHERIC  MILLIMETER  WAVE  PROPAGATION  MODEL 
H.  J.  Liebe* 


The  neutral  atmosphere  is  characterized  for  the  frequency  range 
from  1  to  300  GHz  as  a  nonturbulent  propagation  medium.  Attenuation 
and  propagation  delay  effects  are  predicted  from  meteorological  data 
sets:  pressure,  temperature,  humidity,  suspended  particle  concentra¬ 
tion,  and  rain  rate.  The  physical  data  base  of  the  propagation  model 
consists  of  four  terms:  (a)  resonance  information  for  30  water  vapor 
and  48  oxygen  absorption  lines  in  the  form  of  intensity  coefficients 
and  center  frequency  for  each  line;  (b)  a  composite  (oxygen,  water 
vapor,  and  nitrogen)  continuum  spectrum;  (c)  a  hvdrosol  attenuation 
term  for  haze,  fog,  and  cloud  conditions;  and  (d)  a  rain  attenuation 
model.  Oxygen  lines  extend  into  the  mesosphere,  where  they  behave  in 
a  complicated  manner  due  to  the  Zeeman  effect.  The  geomagnetic  field 
strength  H  is  required  as  an  additional  input  parameter.  Each  0? 
line  splits  Droportionally  with  H  into  numerous  sub-lines,  whicirare 
juxtaposed  to  form  a  Zeeman  pattern  spread  over  a  megahertz  scale. 
Patterns  for  three  main  polarization  cases  are  calculated.  Detailed 
examples  for  model  atmospheres  provide  basic  millimeter  wave  propaga¬ 
tion  information  over  the  height  range  0  to  100  km  of  the  neutral 
atmosphere.  _ 

Key  words:  atmospheric  attenuation;  delay  effects;  millimeter  wave  properties 
of  air;  propagation  model;  oxygen  Zeeman  patterns 


1.  INTRODUCTION 

This  report  is  both  an  update  and  a  recap  of  a  Millimeter-wave  Propagation 
Model  (MPM)  based  on  geophysical  input  parameters.  The  atmospheric  propagation 
medium  is  described  by  measurable  quantities  for  which  spatial  and  temporal  statistics 
are  assumed  to  be  known.  Laboratory  experiments  (Liebe  et  al.,  1977;  Liebe,  T9S4) 
and  analytical  efforts  (Liebe,  1982;  1983)  conducted  at  the  Institute  for  Telecommu¬ 
nication  Sendees  (ITS)  have  played  a  major  role  in  establishing  an  KPN  which 
relates  easy-to-obtain  meteorological  data  (pressure,  temperature,  humidity,  etc.) 
to  dif fieult-to-measure  propagation  factors  (attenuation,  delay,  medium  noise,  etc.) 
in  a  most  direct  manner.  Statistical  problems  have  been  treated  separately  at  ITS; 
for  example,  those  related  to  speeifie  locations  (Allen  et  al.,  1983),  to  rain 
climates  (Dutton  et  al.,  1983),  or  to  rain  clouds  (Allen,  1983). 


Atmospheric  propagation  limitations  dominate  most  consideratons  in  the  advance¬ 
ment  of  millimeter  wave  applications  (Crane,  1981).  Adverse  weather  causes  radio 
signal  degradations  due  to  rain,  wet  snow,  suspended  particles,  and  water  vapor. 

The  MPM  is  structured  to  solve  this  complex  problem  in  an  efficient  manner--comp1 i- 
cated  microphysical  methodology  is  reduced  to  simple  radio  engineering  terms.  Such 
propagation  modeling  provides  a  cost-effective  means  of  predicting  the  performance 
of  a  system  for  its  intended  use.  A  propagation  model  can  take  into  account  all 
performance- influencing  factors  of  the  atmosphere  which  in  the  actual  operating  en¬ 
vironment  would  be  difficult,  if  not  impossible,  to  identify. 

The  propagation  model  is  a  conglomerate  of  individual  routines.  Each  routine 
presents  the  most  practical  simulation  approach  to  a  particular  propagation 
hindrance  and  has  been  validated  by  calibrated  measurements  and/or  theory.  Milli¬ 
meter  wave  signal  degradations  are  described  by  theories  of  molecular  resonances, 
Rayleigh  absorption,  Nie  scattering  loss,  and  Debye  absorption.  One  role  of  the  MPM 
is  to  express  these  theories  or  suitable  approximations  in  terms  of  environmental 
factors,  such  as  pressure,  temperature,  humidity,  suspended  particle  concentration, 
rain  rate,  and  geomagnetic  field  strength. 

The  MPM  basic  to  this  report  was  given  by  Liebe  (1981),  and  a  recap  of  input/ 
output  parameters  is  presented  in  Section  1.1.  The  model  has  found  many  practical 
uses,  as  refleeted  in  Section  1.2  (Allen  et  al.,  1983;  Oamosso  et  ah,  1982; 
Pierluissi  et  al.,  1982;  and  others).  The  heart  of  the  MPM  consists  of  more  than 
400  spectroscopic  parameters.  Some  of  these  have  been  revised  to  take  into  consider¬ 
ation  new  research  findings  (Burch,  1982;  Davies  and  Oli,  1973;  Endo  and  Mirushiaa* 
1982;  Flaud  et  al.,  1981;  Hartmann  and  Kunti,  1983;  Messer  et  al.,  1983;  Pickett 
et  al.,  1981;  Rosenkrana,  1982;  Thomas  and  Nordstrom,  1982).  The  updated  MPM  also 
contains  a  rain  attenuation  term  (Olsen  et  al.,  1978;  Stutsman  and  Oishman,  1982). 
Several  physical  models  of  atmospheric  millimeter  wave  properties  provided  references 
for  the  MPM  reported  here  (Waters,  1976;  Zrashevskiy,  1976;  Falcone  et  al.,  1979; 
Misushiaa,  1982a,  1982b;  Smith,  1981;  Smith,  1982). 

1.1  The  Input/Output  Parameters 

The  Hi  1 1 imeter-wave  Propagation  Model  (MPM)  calculates  radio  path  parameters 
defined  in  this  section  based  on  an  earlier  paper  (liebe,  1981).  For  a  detailed 
description  of  some  of  the  quantities  involved  and  for  a  discussion  of  the  under¬ 
lying  theory,  the  reader  is  referred  to  that  paper.  Germane  to  the  propagation 
model  is  a  maeroscopie  measure  of  the  interaction  between  radiation  and  absorbing 
species  in  the  form  of  the  complex  refractivity  N  (in  parts  per  million). 


ppm 


(1) 


N  =  N0  +  0(f)  +  jN"(f ) 

The  refractivity  consists  of  a  frequency- independent  refractivity  NQ  plus  various 
spectra  of  refractive  dispersion  D(f)  and  absorption  N"(f).  In  radio  engineering 
it  is  customary  to  express  the  imaginary  part  of  (1)  as  specific  power  attenuation 
a  and  the  real  part  as  propagation  delay  0  (with  reference  to  vacuum);  that  is. 


and 


a 

»  0.1820f  N**(f) 

dB/km 

or 

(2) 

a* 

=  a/4.343 

Np/km 

* 

0 

«  0.02096f[NQ  +  0(f)] 

radians/ km 

or 

(3) 

0* 

=  B/2rtf 

ns/ km 

♦ 

where  the  frequency  f  is  in  gigahertz  (GHz)  throughout. 

The  atmosphere  is  characterized  for  the  frequency  range  1  to  300  (1000)  GHz  as 
a  nonturbulent  propagation  medium.  The  medium  is  treated  as  an  ensemble  of  molecules 
and  particles  in  which  spectral  features  of  water  vapor,  oxygen,  suspended  water 
droplets  (hydrosols),  and  precipitation  constitute  the  main  absorbers.  An  adequate 
permanent  spectroscopic  data  base  for  the  calculation  of  N  consists  of  seven  addi¬ 
tive  parts: 


•  the  frequency- independent  refractivity  of  moist  air  given  by  (Liebe,  1981; 
Hill  et  al..  1982) 


«  (2.5S9p  +  4l.6ev  *  2.39e)9  ppm 


(4) 


*  resonance  information  for  30  H^0  lines  between  22  and  997  GHz  (N^), 

•  resonance  information  for  48  0,  lines  between  49  and  834  GHz  (N. ), 

4,  W 

(*)  optional;  resonance  information  for  trace  gas  lines  (Section  3.2)  (Nf>  ^ 

•  continuum  absorption  due  to  >  $QG  H,0  lines  above  1000  GHz  (hi, 

&  W 

•  continuum  due  to  nonresonant  0^  and  pressure- induced  Nj  absorption  (M^), 

*  suspended  water  droplet  (haze,  fog,  cloud)  attenuation  tern 

♦  two-coeff ieient  rain  attenuation  model  (Olsen  et  al.,  1978)  (Nft). 


) 


The  physical  state  of  the  auaosphere,  which  determines  the  absorber  population  in  a 
radio  path,  is  described  by  the  HPM  input  parameters: 
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p 

p 

T 

0 

e 

v 


dry  air  pressure  in  kPa  (1  kPa  =  10  mbar),  where  the  measurable 
barometric  pressure  is 

p  +  e  kPa  ;  ( H '. } 

temperature  in  degrees  Kelvin  ( K) ,  which  for  the  MPM  is 

reduced  to  a  relative  inverse  temperature  parameter 

300/ T  ;  (Fr ) 

water  vapor  partial  pressure,  where  vapor  concentration  is 

7.21 9e0  g/m3  (6a) 

and  both  e  and  v  are  limited  by  relative  humidity 


RH  «  (e/es)100  =  5.750  v0-6  *  IG^9,8346  ’  10>  <  100%  (6b) 

(e  is  the  saturation  pressure)  ; 

5  3 

w  suspended  water  droplet  concentration,  g/m  ;  and 

R  point  rain  rate,  ram/h 


The  dispersion  and  absorption  spectra  are  formulated  from 


o(f)  *  l  (sf*).  ♦  n*  ♦  I  (sf*).  ♦  n;  +  n;*  N* 


i*l 


and 


N“(f )  =  l  (S F“)<  ♦  IIJ  ♦  J  (SF“).  ♦  N“  ♦  N;  ♦ 


i»l 


optional 

n„ 


♦  |  (SF‘)i  ♦  ...  pp®  (7) 


♦  J  (SF-L  «...  ppm  (8) 
<  * 


dry  air  |  humidity  fog, rain  {  trace  gases 


where  $  is  a  line  strength  in  kite  and  F\  F°  are  real  and  imaginary  parts  of  a  line 
shape  function  in  GHt'*.  Both  expressions  are  discussed  in  detail  in  Sections  2 
and  3  where  the  contributing  lines  are  listed  in  Table  1  (n^  *  48  for  0^,). 

Table  2  (nfa  »  30  for  H^O),  and  Table  6  (n£  a  2  for  CO  and  n^  »  32  for  0^). 

For  illustrative  purposes,  it  is  helpful  to  introduce  Figure  1,  which  presents 
an  overview  of  the  clear  air  (w  »  0,  R  8  0)  portion  of  the  NPH  applied  to  predict 
specific  attenuation  a(f)  and  refractive  dispersion  0(f)  for  sea  level  conditions 
at  various  relative  humidities.  Significant  attenuation  extends  over  several  orders 
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of  magnitude  (<10  to  >10  dB/km).  Outside  the  oxygen  resonance  regions  (see 
curve  5),  the  attenuation  rate  is  approximately  proportional  to  absolute  humidity  v 
(6a)  and  is  highly  frequency-dependent.  The  low-frequency  wing  of  the  dispersion 
spectrum  from  the  strong  557  GHz  HjO  line  makes  its  influence  known  down  to  about 
100  GHz.  The  millimeter  wave  range  (30  to  300  GHz)  displays  four  "windows"  with 
relatively  low  values  of  a.  The  frame  for  these  windows  is  set  up  by  resonance 
absorption  due  to  oxygen  (centered  at  60  and  119  GHz)  and  due  to  water  vapor  (centered 
at  22,  183,  and  325  GHz). 

Delay  (3), (7)  and  absorption  (2), (8)  effects  of  the  propagation  medium  are 
functions  of  individual  absorber  abundances  («  p/T,  v,  w,  and  R)  along  a  radio  path, 
and  of  frequency  f.  In  addition,  emission  originates  from  each  absorber  volume 
element  in  order  for  the  medium  to  remain  in  local  thermal  equilibrium.  Radio  path 
measures  for  these  three  interrelated  effects  are: 

a)  total  attenuation  in  units  of  decibels  (dB)  or  nepers  (Np) 


L 

A  =  / a(x)dx  d3  or 

o 

A*  =  A/4.343  Np 


(9) 


b)  total  delay  in  units  of  radians  (rd)  or  nanoseconds  (ns) 
L 

B  3  /B(x)  dx  rd  or 

o 


and 


B*  =  B/2ttf 


ns 


(10) 


c)  sky  brightness  temperature  (Smith,  1982)  in  units  of  degrees  Kelvin  (K) 

**■> 

T g  3  /  ci*(x)T (x)  exp[-A*(x)]dx  ♦  £.9  exp(-A#)  K  ,  (11) 

where  dx  is  an  increment  of  the  path  length  L  in  km,  T ( x )  is  the  kinetic  temperature 
distribution  along  the  path,  and  hQ  is  the  initial  (observer)  height.  Path  attenua¬ 
tion  A  quantifies  the  amount  of  energy  extracted  from  a  plane  wave  propagating 
through  the  atmosphere,  delay  b  is  a  measure  of  the  excess  traveling  time  with  ref¬ 
erence  to  vacuum,  and  brightness  T&  indicates  noise  emission  (between  2.9  K  for 
A*  «  0  and  T  .  .  for  A*  a  -)  from  the  medium  to  the  starting  level  h  . 

dtUD\6f\t  v) 


6 


?J 


is 


?N 


i: 


The  cumulative  radio  path  quantities  (9)  to  (11)  are  the  primary  output  of  the 
MPM.  They  are  evaluated  for  paths  traversing  a  spherically  stratified  (concentric 
layers)  atmospheric  extending  in  height  between  h0  and  hf.  The  final  height  is 
determined  by  the  eventual  absence  of  absorbers.  Typically,  hf  =  100  km  for  dry 
air,  30  km  for  water  vapor,  10  km  for  water  clouds,  and  5  to  0  km  (height  of  0°C 
isotherm)  for  rain.  The  atmosphere  is  structured  in  discrete  height  levels  (e.g., 
layers  between  h  =  0  and  30  km),  each  identified  by  a  set  of  p(h),  T(h),  and  RH(h) 
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data  that  are  used  to  calculate  specific  values  of  a,  3,  and  NQ  for  the  layer 


Cloud  water  profiles  w(h)  can  be  added  to  saturated  (RH  =  100%)  air  conditions,  as 
well  as  constant  rain  rates  R.  The  MPM  can  draw  p-T-RH  height  profiles  from  a 
catalog  of  ten  model  atmospheres,  and  w(h)  from  several  cloud  models  (Section  4.1). 
"Ray  tracing"  is  used  to  evaluate  the  path  increments  dx.  A  ray  is  assumed  to 


start  from  the  initial  level  hQ  with  an  elevation  angle  (measured  from  the  hori¬ 


zontal)  and  proceed  through  the  atmosphere,  gaining  the  height  interval  dh  =  h  -  hQ 
while  being  subjected  to  refractive  bending.  The  traveled  path  differential,  includ¬ 
ing  refractive  bending,  is  given  by 


,  r  /  10°  +  N ‘ ( h  )  \  /  rc  +  h  \  ,2  1 

dx  =  dh/{'-  [  (w )  coH  } 


+  h. 


,2  ,1/2 


km 


02) 


where  r^  =  6357  km  is  a  standard  (45°N)  earth  radius.  Equation  (12)  is  the 
spherical  form  of  Snell ‘s  law  and  is  strictly  valid  only  when  N  is  purely  real,  but 
for  ehe  case  at  hand  it  is  a  good  approximation,  since  the  values  for  N"  are  gener¬ 


ally  small  compared  with  N'  3  NQ  +  D. 


The  shortest  way  through  the  atmosphere  is  the  zenith  path  i>  =  90°,  where 
dx  =  dh.  A  numerical  integration  through  the  layered  atmosphere  implies  simply  a 
summation  of  the  specific  contributions  at  each  1c. cl,  applying  Simpson's  rule. 


Clear  air  zenith  attenuation  A2  is  said  to  be  causeo  uy  "one  air  ..ass."  To  first 


order,  the  air  mass  for  slanted  radio  paths  down  to  tj>  >  10°  increases  according 
to  the  secant  law,  dx  a  dh/sin  ty,  which  has  been  experimentally  verified  for  slant 
path  attenuation  A$  (Altshuler  et  al.,  1978).  The  behavior  of  slanted  (^  >  10°) 
radio  paths  can  be  predicted  from 


As  =  Aj/itn  y 


and 


3$  3  B^/sin  1 1> 


(13) 


« 

* 


9i 


< 

> 


? 
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Only  low  elevation  angles  (i|/  <  10°)  require  a  detailed  evaluation  of  (12). 
Refraction  is  most  pronounced  in  the  lowest  layers  of  the  atmosphere,  which  have  to 
be  divided  into  many  sublayers  for  an  accurate  tracing  of  the  ray  path  (Gallery 
et  al.,  1983).  For  ip  <  1°  and  values  of  h  close  to  h  ,  Equation  (12)  becomes 
numerically  unstable  and  approximations  of  the  two  bracketed  terms  have  to  be  made 
by  series  expansions.  At  this  point,  path  increments  dx  become  very  sensitive  to  a 
particular  height  distribution  of  absorbers.  For  example,  the  tangential  (<p  =  0°) 
air  mass  through  the  U.  S.  Standard  Atmosphere  (NOAA,  1976)  for  dry  air  is  38  times 
the  zenith  value  (Gallery  et  al.,  1983),  while  for  water  vapor  the  tangential  mass 
is  about  70  to  180  times  the  zenith  value. 

Three  useful  cumulative  measures  describing  the  propagation  medium  are: 


a)  the  refractive  path  delay 


Bo  /  N0(x)dx 

o 


mm 


or 


B*  =  3  336  Ba 
o  o 


PS 


(14) 


b)  the  total,  nath-integrated  water  vapor 


V  =  /  v(x)dx 
o 


mm 


and 


05) 


c)  the  total  suspended  dropiet  water 


I 

I 


/  w(x)dx 
o 


mm 


(16) 


The  MPM  operates  as  follows: 


INPUT 


S 


frequency  range  f^>_»  GHz)  to  ?2(  300  GH2), 

which  may  be  extended  to  1000  GHz,  although 
purely  on  theoretical  grounds 

horizontal  path  of  length  L  or  slant  path  with 
elevation  angle  ">  at  hQ 

pressure  and  temperature  p,fl 

water  vapor  concentration  v(RH) 

hydrosol  concentration  w 

rainfall  rate  R 


(Fig.  1) 


(12), (13) 
(5a, 5b) 
(6a, 6b) 


$ 

w 


OUTPUT  specific  rates  a(f),8(f)  (2), (3)  , 

cumulative  measures  A(f),B(f)  (9), (10) 

medium  noise  Tg(f)  (11) 

medium  descriptors  BQ,V,W  (14). (16) 


The  heart  of  the  model  is  a  complex  refractivity  N  (7), (8)  expressed  in  measurable 


quantities.  The  body  of  this  report  consists  (a)  of  improvements  and  additions  to 
the  formulations  for  N  over  two  separate  height  ranges,  h  =  0  to  30  km  (Section  2) 
and  30  to  100  km  (Section  3);  (b)  of  a  discussion  of  the  MPM  computer  routine 
(Section  4);  and  (c)  of  samples  of  representative  results  in  numerical  and  graphical 
form  (Section  5  and  Appendix  A). 

1.2  Model  Applications 

Since  the  introduction  of  the  first  MPM  version  (Liebe,  1981),  the  model  has 
proven  useful  for  a  variety  of  applications  in  millimeter  wave  communications, 
remote  sensing,  and  radio  astronomy.  Tape  copies  of  the  model  together  with  program 
listings  for  PI  and  P2  (see  Section  4)  have  been  requested  by  51  users  from  industry 
(23),  Government  (13),  universities  (7),  and  foreign  countries  (8).  The  Institute 
(ITS)  has  worked  with  the  MPM  to  lend  support  to  the  following  applications,  to  name 
but  a  few: 

•  to  develop  simple  models  for  predicting  transparency  in  the  four  atmospheric 
millimeter  wave  ranges  (Liebe,  1983); 

•  to  present  selected  examples  of  the  unique  transfer  characteristics  of  dry  air 
in  the  50  to  70  GHz  frequency  range  (Liebe,  1982); 

•  to  estimate  location-dependent  attenuation  between  10  and  100  GHz  for  major 
U.  S.  population  centers  (Allen  et  al.,  1983); 

•  to  plan  the  development  of  a  system  performance  model  for  terrestrial  milli¬ 
meter  wave  communications  (Liebe  et  al.,  NTIA  Technical  Memorandum  83-93, 

August  1983,  limited  distribution); 

•  to  provide  numerical  data  on  the  effective  height  concept  h$  =  A/a  (Section 
5.2)  over  the  range  5  to  400  GHz  using  three  model  atmospheres  (U.  S.  Standard, 
Arctic  Winter,  and  Tropical  Summer)  for  CCIR  Study  Group  5,  International 
Working  Party  3,  July  1983  (F.  Fedi,  Chairman); 

•  to  evaluate  water-vapor-dependent  attenuation  data  from  laboratory  and  field 
exper  iments  conducted  at  20.6,  31.6,  96.1,  110,  and  337  GHz  as  suimiarized  in 
Appendix  C; 


•  to  compare  the  propagation  characteristics  over  the  range  30  to  150  GHz  for 

two  sites  selected  at  h  =  1  km  in  Japan  and  at  h  =  3  km  in  Spain  for  millimeter 
wave  radio  astronomy  (M.  Lebenbaum,  private  communication,  1983); 

•  to  interpret  data  from  a  radiometric  observation  of  the  25"  0^  line  (see 
Section  3.1  and  Figure  A2)  performed  over  Europe  from  an  aircraft  flying  at  a 
height  of  10  km  (Hartmann  and  Kunzi,  1983). 

2.  ABSORPTION  IN  THE  LOWER  ATMOSPHERE  (0  to  30  km) 

Strong  interactions  do  occur  between  millimeter  waves  and  moist  air  as  evident 
from  Figure  1.  Other  phase  states  of  H2O  also  act  as  strong  absorbers  in  a  radio 
path,  each  with  a  behavior  that  is  distinctly  different  across  the  millimeter  wave 
spectrum.  Specific  attenuation  for  five  H~0  states  is  illustrated  in  Figure  2, 

l  22 

assuming  approximately  1  mm  absorber  thickness  («3.4  *  10  H20  molecules).  Param¬ 

eterization  of  the  theories  that  need  to  be  applied  to  express  the  individual  spectra 
of  (8)  in  terms  describing  atmospheric  conditions  is  outlined  in  this  section  includ¬ 
ing  some  revisions  with  regard  to  the  initial  model. 

2.1  Molecular  Effects 

Gaseous  02  and  H20  are  considered  to  bw  the  principal  absorbers  in  moist  air. 
Pressure  broadening  of  spectral  lines  leads  to  two  different  types  of  frequency 
responses;  namely,  sharp  resonance  lines  and  continuum  spectra,  which  are  treated 
separately  in  the  following  for  atmospheric  conditions  up  to  heights,  h  =  30  km. 

Self-  and  foreign-gas-broadening  influences  have  to  be  taken  into  account. 

2.1.1  Local  Line  Absorption  and  Dispersion 

For  the  shape  of  an  absorption  line,  the  Van  Vleck-Weisskopf  function  (Kemp, 
1979)  as  modified  by  Rosenkranz  (1975)  was  chosen  in  the  form 

f  /Y  -  (V  -  f)<5  Y  -  (v  +  f)6  .  , 

F*‘  =  —  ( - 2 -  + - - 2 - ]  GHz  .  (17) 

v0  \(v0  -  f)2  +  Y2  (v0  +  f)2  -  Y2  / 

A  Hilbert  transform  of  (17)  yields  the  cot  responding  dispersion  profile  (D.  Brown, 
private  communication,  1982;  R.  Hill,  private  communication,  1982) 
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CURVE  I  kCDIUM 


WATER,  Imm 
RAIN,  20  mm/h 
CLOUO  WATER 
MOIST  AIR,  101  kPo 
WATER  \AP0R 


I  g/m1  i  I  mm/ km 


300°  K 


BJ \±L 

w  ill L  _ 


SI 


FREQUENCY,  GHz 


*)  Note: 


The  phase  states  (absorption  theories)  are:  1)  water  (Debye 
relaxation),  2)  rain  (Mie  scattering),  3)  hydrosol  (Rayleigh 
absorption,  4)  moist  air  (pressure-broadened  spectral  line, 
numbers  shown  indicate  selected  peak  attenuations),  and  5)  pure 
water  vapor  (self-broadered  lines). 


Figure  2.  Specific  HpO  attenuation  over  a  frequency  range  from  5  to 
600  GHz  for  five  phase  states  normalized  to  equal  absorber 
amount  (1  mm). 


f(vQ  -  f)  +  y(y  +  f6)/vQ 
^  (v0  -  f)2  +  Y2 


(v  +  f)  +  Y(Y  -  f6)/v  ?  \  i 

— - -  -  ~  }  GHz"  ,  (18) 

(vQ  +  f)2  +  y2  0/ 


which  is  slightly  different  from  F'  as  given  by  Liebe  (1981). 

The  line  parameters  are  calculated  according  to  the  scheme  below: 


02  1 ines  in  air 

H20  lines  in  air 

Equation 

ID 

a 

b 

(7), (8) 

vQ  ,  GHz 

Table  1 

Table  2 

S  ,  kHz 

a-jpe3  exp[a2(l  -  0)]1O~6* 

b^G3,5  exp[b2(l  -  0)] 

09) 

Y  »  GHz 

(0.8-a,)  3 

a3(pe  0  +  1.  Ie8)10~'3 

b3(p00,8  ♦  4.80e6)10"3 

(20) 

6  ,  1 

a5 

a4p0 

0 

(21) 

* 

See  Appendix  B 


Line  center  frequencies  v  and  the  spectroscopic  coefficients  a^  to  a^,  b^  to 
for  strength  S  and  width  y  are  listed  in  Tables  1  and  2.  Revisions  have  been 
made  to: 


vo(02) 

41 

i3(e) 

a6 

W 


(Endo  and  Mizushima,  1982)  , 

(Appendix  B)  , 

(Set zer  and  Pickett,  1977)  , 

(Pickett  et  al.,  1981)  , 

(Messer  et  al . ,  1983  )  , 

(Flaud  et  al.,  1 981;  Mizushima,  1982a)  , 


b 


3 


The  oxygen  coefficients  to  are  based  largely  on  solid  experimental 
evidence  (Liebe  et  al . ,  1977).  The  spectroscopic  fundamentals  for  the  molecular 
calculation  procedure  of  a^  and  a^  are  outlined  briefly  in  Appendix  B.  Even  after 
drastic  simplifications,  there  remain  336  spectroscopic  coefficients  (Table  1)  to 
describe  for  an  MPM  the  detailed  line  spectrum  of  dry  air.  The  temperature  depend¬ 
ence  of  the  overlap  coefficients  a4  evolves  from  a  complicated  iterative  procedure 
(Rosenkranz,  1975)  that  was  evaluated  over  the  temperature  range  T  =  200  to  300  K. 
The  results  were  approximated  by  a^  (21).  The  even  simpler  approximation  a^  *  2 
was  attempted,  and  a  comparison  of  both  with  rigorous  calculations,  as  shown  for  an 
example  in  Table  3,  resulted  in  opting  for  a5  as  listed  in  Table  1.  A  more  rigorous 
theory  for  a^.a^  has  been  developed  by  Smith  (1981)  based  on  our  data. 

The  spectroscopic  water  vapor  coefficients  b^  also  have  a  molecular  base  in 
Equation  (B2).  The  asymmetric  rotor  molecule  possesses  a  strong  electric  dipole 
moment  (u)  and  requires  different  assumptions  about  selection  rules  (vQ),  matrix 
elements  (|<J>|),  energy  levels  (E£),  partition  function  (Q),  and  nuclear  spin  (g^). 
The  line  strength  coefficients  b^  were  compared  with  several  reported  H^O  line  data 
bases  (Rothman  et  al.,  1983a;  Poynter  and  Pickett,  1981;  Mizushima,  1982a).  These 
data  bases  differ  slightly  from  one  another  and  the  b1  values  given  in  Table  2  are 
derived  from  Flaud  et  al.  (1981).  linewidth  data  b3  for  the  22  GHz  and  183  GHz 
lines  are  based  on  controlled  laboratory  experiments;  the  remaining  ones  were 
approximated  following  the  theoretical  work  by  Davies  and  Oli  (1978).  The  MPM 
version  of  the  atmospheric  water  vapor  line  spectrum  requires  120  spectroscopic 
coefficients  (Table  2)  since  overlap  corrections  need  not  be  made.  Water  vapor 
lines  are  reported  to  display  small,  vapor-pressure-induced  shifts  of  tl.'  center 
frequencies  vQ,  which  are  neglected. 

2.1.2  Continuum  Absorption 

The  continuum  spectrum  in  (7)  and  (8)  consists  of  dry  air  and  water  vapor 
absorption  terms, 

N"  »  N"  ♦  N“  ,  (22) 

c  p  t? 

while  the  respective  refractive  dispersions  are  neglected;  i.e.,  ■-  0  and  0. 

The  continuum  varies  monotonically  with  frequency  and  needs  to  be  added  to  the 
selected  group  of  local  and  H^O  resonance  lines  (Tables  1  and  2)  in  order  to 
correctly  predict  atmospheric  millimeter  wave  attenuation  in  the  window  ranges 
in  between  lines. 
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Table  1.  Revised  Oxygen  Line  Parameters 
(compare  with  Liebe,  1981:  Table  4a) 


££ 

&  o 

al 

a2 

a3 

a4 

a5 

a6 

S  GHz  (kHz/kPa)10‘ 

6 

MHz/kPa 

1/kPa 

ft'  49. 4523  T9 

0.12  - 

-JX.A3C 

d.40 

5.60  _ 

1*7 

a 

ft  49.9622  57 

0.34 

10. 

720 

8.50 

5.60 

1.7 

0 

S';  50.474238 

0.94 

9, 

690 

8.60 

5.60 

1.7 

0 

50.967743. 

2*46 

.  8. 

690 

8. 70 

5.50 

1.7 

_ a 

■  51.503350 

6.08 

7. 

740 

8.40 

5.60 

1.8 

0 

K<  52.0214C9 

14.14 

6. 

*40 

9.23 

5.50 

j.e 

0 

W  52.54.2393 

31.02  . 

6, 

000 

9.40  . 

.  .5*.?.Q _ 

.1,8 

..  0 

<\>  53.066906 

64.10 

5. 

220 

9.70 

5.30 

1.9 

0 

fc  53.595746 

124.70 

4, 

4  90 

1C. CO 

5.40 

1.8 

0 

&  54,129999 

_  228.00 

3, 

9 10 

10.20 

4,8.0. 

_2*Q 

a 

1  54.671157 

391,80 

3. 

190 

10.50 

4.83 

1.9 

0 

Co  55.221365 

631.60 

2, 

a20 

10.79 

4.17 

2.1 

0 

ft  55.7836C0 

953.50 

2, 

115 

U.*0 

3.75 

. 2.1 

0 

sO  56.264777 

548.90 

0, 

010 

16.46 

7.74 

0.9 

3 

ft  56.363387 

1344.00 

1. 

655 

11.44 

2.97 

2.3 

0 

£  .  56*968180 

1763.00 

1, 

251 

11.81 

....  2,12 _ 

-2*5 

„.Ql 

9  57.612481 

2141.00 

0. 

910 

12.21 

0.94 

3.7 

0 

&  56.323874 

2386.00 

0. 

621 

12.66 

-0.55 

-3.1 

0 

SJ  58.446589 

1457.00 

0. 

3  79 

14,49 

5.97 

0.8 

0 

V  59.164204 

2 404.00 

0. 

386 

13.19 

-2.44 

0,1 

0 

&  59.590982 

2112.00 

0. 

>07 

13.60 

3.44 

0.5 

0 

^  60, 3060 57 

2124,00 

0, 

207 

13,82 

-4,13* 

Q.7 

Q 

1  60.434775 

2461,00 

0. 

386 

12.97 

1.32 

-1.0 

0 

61.150558 

2504.00 

c. 

621 

12.48 

-0.36 

5.8 

0 

to  61.800152 

2208.00 

0. 

910 

12.07 

-1.59 

2.9 

0 

62.411212 

1933.00 

1. 

255 

11.71 

-2.66 

2.3 

0 

£\  62.486253 

1517.00 

0, 

078 

14,68 

-4.77* 

0,9 

0 

£  62.997974 

1503.00 

1, 

660 

11.39 

-3.34 

2*2 

0 

m  63.563515 

1037.00 

2. 

110 

11.08 

-4,17 

2.0 

0 

R  64.127764 

733.50 

2, 

620 

1C. 78 

-4.48 

2.0 

0 

E'\  64,678920 

463.53 

3. 

190 

10.50 

-5.10 

1*8 

0 

p?  65,224067 

274.80 

3* 

310 

10.20 

-5.10 

1.9 

3 

K>  65.764769 

153.00 

4, 

430 

10.03 

-5.70 

1.8 

0 

66.302083 

80.09 

5. 

229 

9.70 

-5,50 

1.8 

0 

M  66.836827 

39,46 

6  1 

000 

9.40 

-5.90 

1.7 

0 

[*:•;  67.369595 

IS,  32 

6, 

140 

9,20 

-5.60 

1.8 

0 

■  •>  67.900862 

8.01 

7. 

740 

9,90 

-5.50 

1.7 

0. 

68.431501 

3*30 

8. 

690 

P.70 

-5.70 

1.7 

0 

Bft  68.960306 

1.28 

9, 

690 

8*60 

-5.60 

1.7 

0 

■  /  69*489021 

0.47 

10. 

720 

9.50 

-5.60 

1.7 

0 

Ef  70.017342 

0.16 

11, 

330 

9,40 

-5,60 

1.7 

0 

IS  JL2Ji.U03L4j. 

945.00 

9. 

000 

-JLL 

_ 0*9 — 

0 

D=doublet 


Table  2.  Revised  Water  Vapor  Line  Parameters 
(compare  with  Liebe,  1981:  Table  4b) 


22*235080 

67*813000* 

120*127000* 

183.310117 

321.225644 

325.152919 

336*187000* 

380.197372 


kHz/kPa 

0.1120  2.143 

0.0011  8.730 

0.0007  8.350 

2.4070  0.653 

0.0450  6.160 

1.5910  1.520 

0.0009  9.800 


12.6400 


437.346667  0.0625  5.020 

439.150812  0.9250  3.560 

443*018295  0.1900  5.020 

448,001075  10.9900  1.370 

470.888947  0.3310  3.570 

474.689127  1.3100  2.340 

488. 49H  33 _ L UUQ 

503.568532  0.0354  6.690 

504.482692  0.0118  6.690 


1.020 

7.330 

5.020 

3.560 

5.020 

1.370 

3.570 


620.700807 

658.015000* 


5.3960 

0.4600 


2.340 

7.760 


MHz/kPa 

28.10 

25.00 

_2_7*P0_ 

28.20 

25.50 
_29.0Q_ 

24.00 

28.50 
25.00 

25.50 

25.80 

_j25.50 

27.60 

26.80 

26.50 
26.00 _ 
24.40 
24.40 

_3U7<1. 

27.00 

32.80 


841.073593 

859.865000* 


0.0123 

0.0150 


8.110 

7.990 


25.40 

29,00 


902.555000* 
906.205524 
916.171582 
970.315022 
V87. 926764 


0.0640 

0.1820 

8.9180 

9.6880 

145.6000 


8.350 

5.040 

jum. 

1.840 

0.180 


27.80 

25.40 

27.90 
26.60 

29.90 


*  First  vibrational l,y  excited  line  (Flaud  et  al.»  1981) 
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Table  3. 


f>f -?U£.'Ct 

<  uH’l 


♦  0.0 
4  2.5 

♦  5.<; 

♦  ».? 

70.0 

72. 7 
57. C 

77.7 

60.0 

72.7 
<7.0 
7  7.7 
70.0 

72.7 
77.0 

77.7 
80.0 

82.7 
87.0 

87.7 
70.0 
*>7.5 
q7 ,0 
9  7,7 

100.0 

102.7 
107.0 
10  7.7 
mi.L 
112.? 
117,0 

117.7 

120.0 

122.7 
127.0 
12  7. S 
UO.C 

112.7 

1  J7.0 
1  1  7.7 

m.o 


n?7d,fPL0Xiftl0nS,:°r.the  terilPeratur^  dependence  of  Rosenkra 
(1975)  overlap  coefficients  for  the  02  microwave  spectrum. 


1)  The  "a,-  =  2"  Approximation: 


(1ST  A 1  8  AT  270  K,  p:  101  kPA 


AT  T  £  NUA T ION 
<  08  /*  “  ) 


PESCrr.f  fR*C6 
or  a p v  s o * . 


Pp‘7t  riiPfASION 
«  AAO/XK) 


PtRCtNT  f  K  R  OR 

or  appro*. 


.on 

.1  OS 

.1*6 

.*51 

1.0*3 

♦  ,■*•1 
14.349 
20.674 
15.707 
♦  .♦<7 
1.129 
.711 
.110 
.2  !♦ 
» l  7  9 
.12  7 
.102 
.08  8 
.077 

•  0*7 
.071 
.0*8 
.0-7 

•  0  5  7 
.070 
.0-  S 
.01  j 
.117 
.1.7 
.♦13 

1.110 
1  .294 
.ITS 
.1'  1 
.0 ‘>U 
.3*9 
.6*  » 
.01  1 
.0:  7 
.a:  i 


-.04 
-.09 
-  .  1  ♦ 
-.18 
-.17 
.10 
.57 
.01 
-.♦8 
.1.8 
.06 
-.90 
-1.28 
-1.39 
-1,44 
-1.45 
-1.41 
-1.33 
-1  .20 
-1.01 
-.77 
-.42 
-.01 
.48 
1  .02 
1 .79 
;  .ii 
2.71 
•'  « 1  4 
2  .44 
1.78 
.47 
-.8? 

-  .  .  76 
-4.40 
-f  .21 
-7.94 
.56 

-11.01 
-l.  .32 
-13.76 


.115 
.172 
•  2  0  / 
.291 
.♦27 
.fcM 
1.177 
1.01  1 
-.173 
-1.6C7 
-1.673 
-1.175 
-.«;  t 
-.794 
-.713 
-.6  S 
-.*2  0 
-.7  9* 
-.3  72 
-.‘36 
-.54* 

-.733 

-.328 

-.5.-2 
-.‘•18 
-  . :  13 

IS 
-.304 
-.4  96 

-.464 

-.460 

-.474 

-.474 

-.644 

-.m 
-.*1 P 

-.46* 
-.if  ‘ 

•  •  f  0  6 
-.6.9 
-.61< 


2)  The  “a5"  Approximation  (Table  1,  fq.  21): 

♦!*A  •??  .in 

*  *81  .»%? 

•1  .6?  ,207 


.6  <*  1 
.los 
.1*6 
.8*2 
.*31 

1.  i“4 

4,*H 

20.-.  74 
17.7b* 


*<■  .  •, 

,r? 

■>*  • : 

,3r.6 

ft?.-' 

.  5  /  3 

S3  .  . 

.<  *  7 

♦2.3 

.4-  1 

S3,  a 

«v8s6 

S7.3 

.654 

tit  .5 

.08  7 

i  .*7.8 

«Jt  t 

1  ?  '  •  5 

. :  *  »■ 

U7.‘ 

.**>  J 

11-’. 1 

.115 

it’.*. 

,1-9 

117.0 

.4  1  i 

tt/.‘ 

1  .  ‘  5  $ 

l?8.c 

t..*4 

1  2  7  *  * 

,1  <9 

l  7  *• .  6 

.H  1 

12».3 

iu.e 

.♦27 
.691 
1.1*7 
1.01  ) 
-.177 
-1,607 
-1.653 
—  1.17* 
-.♦?< 
794 
-.713 
•«6‘  7 
-.620 
-.393 
-.372 
-.534 
-.344 
-.333 
-.3  24 
-.72? 
-.‘if 
-.‘il 
-.863 

-.364 
-.49* 
-.48  4 
-.46  S 
-.434 
-.b*  * 
-.644 
-.b?l 

-.no 

•  *  t  0< 

-.665 

-.66* 

*.6QS 
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The  dry  air  continuum, 

Np  "  Ca0pe2y0/(f2  +  y2}  +  anp202,5]f  ppm  , 


(23) 


makes  a  small  contribution  at  ground  level  pressures  due  to  the  r.onresonant  02 
spectrum  below  10  GHz  (Rosenkranz,  1975)  and  the  pressure- induced  N2  spectrum 
effective  above  100  GHz  (Oagg  et  al . ,  1982).  The  scale  factors  in  (23)  are 


=  6.14  x  io-4 

ppm(kPa  GHz)-1  , 

=  1.5  x  10'10 

ppm(kPa)"2  (GHz)'1 

(23a) 


i 


and  the  width  parameter  (see  Eg.  20), 

=  5.6  x  10'3  (p  +  l.le)0°‘8  GHz/kPa 
o 


(23b) 


for  tne  Debye  spectrum  of  Q?,  is  chosen  to  agree  with  overlap  measurements 
(Liebe  et  al.,  1977)  following  a  suggestion  by  Rosenkranz  (1982). 

The  water  vapor  continuum, 

NJJ  8  [bfepe2,5  ♦  b5e2a3*5]f  ppm  ,  (24) 


has  been  a  major  source  of  uncertainty  in  predicting  millimeter  wave  attenuation 
rates,  especially  in  the  four  window  ranges  (Waters,  1976;  Rice  and  Ade,  1979; 

Emery  et  al.,  1980;  Crane,  1981;  Zammit  and  Ade,  1981;  Burch,  1982;  Zammit  et  al., 
1982;  Liebe,  1983).  Recent  laboratory  experiments,  employing  a  special  high-humidity 
spectrometer  at  f  *  138  GHz,  RH  =  80  to  1Q0S,  and  p  8  Q  to  150  kPa  (nitrogen), 
produced  the  following  results  (Liebe,  1984),  which  are  refleeted  in  the  formulation 
of  (24): 


•  the  scale  factors  in  (24)  for  moist  air  are 

bf  8  1.40  x  io-6  ppm(kPa)"2  GHz"1 
and 

b$  8  5.41  x  io-5  ppm(kPa)"2  GHz"1 


(24a) 

(24b) 
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up  to  saturation  (RH  =  100%),  there  was  no  evidence  of  anomalous  absorption 
at  the  two  test  temperatures,  282  and  300K; 

a  strong  self-broadening  component  b$e2  was  measured  and  tentatively  identified 
to  be  related  to  monomer  H20  line  absorption,  since  its  temperature  dependence 

is  0  •  -  ,  which  is  a  little  higher  than  expected  from  H20  line  broadening 
theory. 

Predictions  with  the  MPM  including  (24)  show  good  agreement  with  absolute  laboratory 
and  field  results  for  specific  attenuation  caused  by  water  vapor  (see  Appendix  C). 

2.2  Suspended  Particle  Effects 

Hydrosol s  in  haze,  fog,  or  clouds  are  the  dominant  suspended  particle  absorbers 
in  the  lower  atmosphere.  Their  refractivity  contribution  is  labeled  Ny  *  (N*  ♦  jir). 
Below  300  GHz,  there  is  relatively  little  effect  on  path  delay  since  *  * 

N;*1.4w  ppm  (25) 

It  is  the  Rayleigh  absorption  approxifftation 

N;  •  V°-,82f  *  «•»  "  «"/[(£•♦  2)2  ♦  (<:•)*]  ppm  (26) 

of  Mie  scattering  losses  (appropriate  for  size  ranges  of  radii  <  50  urn),  which  makes 
a  frequency  dependent  contribution  (Falcone  et  al.,  1979).  TmTdieleetric  data 
e  .  c  of  bulk  water  are  calculated  with  the  Debye  model  reported  by  ehang  and 
Wilheit  (1979),  which  is  valid  for  f  <  300  GHz: 


where 


s  4.9  ♦  0  65  -  113/u )/(!  ♦  (ft)2] 
s  08S  -  n3/o)ft/Cl  ♦  (fO2)  . 


5  9  4.17  *  10°  0  exp(7. 13o) 


Equations  (2),  (26),  and  (27)  are  applied  to  evodel  specific  attenuation  by  hydrosol s 
(haze,  fog,  elouds)  if  their  average  mass  concentration  w(g/ta^)  and  relative 
temperature  o  within  a  radio  path  are  known.  Cloud  ana  fog  conditions  (w  -  0.01  g/m^) 
exist  only  in  nearly  saturated  air  (Rh  *  lOOx);  for  haze  conditions  (w  -  0.01  g/ta"*}. 
the  relative  humidity  taay  be  in  the  90  to  99  percent  range.  Table  4  lists 


numerical  examples  of  o  (dB/km)  for  two  temperatures  (0°  C  and  25°  C)  at  selected 

w 

frequencies. 

Table  4.  Hydrosol  (haze,  fog,  cloud)  attenuation  a 

3  w 

for  a  mass  concentration  w  =  1  g/m 


TEMPERATURE  FREQUENCY  f  ,  GHz 


T 

0 

1 

10 

30 

100 

200 

300 

°C 

dB/km 

0 

1 . 0983 

.001 

.097 

.82 

5.4 

9.3 

10.7 

25 

1 . 0062 

.001 

.051 

.45 

4.2 

10.8 

15.3 

Models  for  c'juds,  fog,  and  haze  are  available  that  describe  spatial  distribu¬ 
tions  of  w(h)  (e.g.,  Chang  and  Wilheit,  1979;  Falcone  et  al.,  1979).  Nonprecipi- 
tating  clouds  and  clouds  producing  stratiform  rain  have  limited  vertical  but,  on 
the  average,  large  horizontal  extents.  Cloud  top/bottom  heights  typically  are  in 
ranges  from  (2  to  8)/(l  to  6)  km,  and  the  liquid  water  concentration  w  may  vary 
between  0.01  and  1  g/m  .  Chemical  depositions  in  clouds  (e.g.,  acid  content)  may 
alter  the  dielectric  properties  (27)  of  pure  water.  Also,  temperature  (roughly 
between  -20  and  +5°  C)  has  some  influence  upon  the  attenuation  (26).  Integrated 
zenith  column  densities  W  (16)  amount  to  values  between  0.05  and  2  mm,  distributed 
in  air  masses  with  water  vapor  column  densities  V  (15)  between  5  and  60  mm.  For 
slant  path  calculations,  both  columns  are  expected  to  display  a  secant  dependence 
(13)  with  elevation  angles  t|>  >  10°. 

The  structure  of  cumulus-type  (convective)  clouds,  such  as  occur  in  thunder¬ 
storms,  is  much  more  complicated.  Here,  distinct  regions  for  updrafts  (large  W) 
and  downdrafts  (large  R),  as  well  as  an  irregular  horizontal  structure,  have  to  be 
distinguished.  An  attempt  to  model  their  influence  upon  millimeter  wave  propagation 
was  reported  by  Allen  (1983).  He  suggests  for  both  types  of  rain-producing  clouds 
simple  correlations  between  hydrosol  concentration  w  and  rain  rate  R  in  the  mean 
(stratus:  w  =  G.7R^,  cumulus-type:  w  -  3  R^‘J)  and  proceeds  to  derive  rain 

cloud  properties  from  readily  available  rain  rate  statistics. 


2.3  Precipitation  Effects 

Precipitation  originates  as  a  hig  ,ly  statistical  event  within  clouds  suspended 
in  saturated  air.  Its  vertical  distribution  is  separated  into  two  regions  by  the 
height  hy  of  the  C°  C  isotherm,  which  can  vary  between  6  km  and  ground  level 
depending  on  season  and  latitude.  The  lower  part  is  mostly  liquid  drops,  and  the 
upper  region  consists  of  frozen  ©articles  with  occasional  supercooled  droplet-load¬ 
ings  by  strong  updrafts.  Widespread  steady  rain  occurs  more  uniformly  and  favors 
small  drop  sizes  (<_  1  mm  diameter)  which  stay  in  the  air  longer.  Heavy  showers  are 
more  localized,  favor  larger  drops,  and  occur  less  frequently. 

The  refractivity  spectrum  of  rain  is  identified  in  (7)  and  (8)  by 
Nr  -  +  JNr  .  Drop  diameters  (0.1  -  6  mm)  and  millimeter  wavelengths  are  compa¬ 

rable,  thus  causing  appreciable  interactions  due  to  Mie  absorption  and  scattering. 

A  very  crude  approximation  for  delay  effects  yields 


N,)  as  0.06  R  ppm  for  f  <  10  GHz  and 

(28) 

Nr  «  0.6R/f  ppm  for  f  >  10  GHz  . 

The  absorption  spectrum  is  formally  expressed  by 


NJ«aR/0.182f  ppm  . 

Bypassing  elaborate,  lengthy  Mie  calculations  (e.g.,  Falcone  et  al.,  1979),  which 
require  drop  shape  and  size  distributions,  as  well  as  the  complex  dielectric 
constant  of  water  (27),  Olsen  et  al.  (1978)  approximated  the  specific  rain  attenu¬ 
ation  by  the  power  law  model 

4r  *  u(f)Rv*fJ  dB/km  .  (29) 

The  frequeney-dependent  coefficient  u  and  exponent  v  were  calculated  using  dropsize 
spectra  of  Laws  and  Parsons  and  a  temperature  of  T  =  0°  C.  A  regression  fit  to 
individual  (u.v)-pairs  over  the  frequency  range  from  1  to  1000  GHz  resulted  in  the 
following  calculation  scheme  (Olsen  et  al.,  1978): 


1  to  2.9  6.39  x  10 

2.9  to  54  4.21  x  10 

54  to  180  4.09  x  10 


si 


V  = 


c  4 


(29a) 


-5 

-5 

-2 


180  to  1000 


3.38 


2.03 

2.42 

0.699 

-0.151 


GHz 

1  to  8.5 
8.5  to  25 
25  to  164 
164  to  1000 


0.851 

1.41 

2.63 

0.616 


0.158 

-0.0779 

-0.272 

0.0126 


An  evaluation  of  (29)  for  point  rainfall  rates  R  between  0.25  (light  drizzle)  and 
150  (severe  thundershower)  mm/h  is  shown  in  Figure  3.  The  validity  of  this  approach, 
which  reduces  theoretical  complexity  to  empirical  simplicity,  is  fairly  well 
established  (Stutzman  and  Dishman,  1982). 

Rain  can  be  a  serious  attenuator  at  frequencies  above  20  GHz;  however,  experi¬ 
ence  is:  "it  rains  only  occasionally."  Actually,  events  of  serious  consequences 
(aR>l  dB/kn()  occur  typically  less  than  one  percent  over  the  period  of  a  year. 

Long-term  statistics  of  rain  attenuation,  however,  are  beyond  the  scope  intended 
for  the  MPM.  Statistical  estimation  methods  and  models  for  the  occurrences  of  rain 
attenuation  are  discussed  in  detail,  for  example,  by  Crane  (1980,1981)  and  by 
Dutton  et  al.  (1983).  Rain  cell  extent  LR  and  path-averaged  rain  rate  R  are  obtained 
from  cumulative  distributions  of  path  attenuation  A  or  point  rate  R  at  many  locations 
over  long  time  periods,  preferably  several  years. 

The  MPM  is  capable  of  treating  rain  identified  by  a  point  rate  R  (mm/h)  meas¬ 
ured  somewhere  along  the  path  extending  over  a  horizontal  length  L  (km).  Then  an 
average  rate  R,  which  is  calculated  with  an  empirical  model  proposed  by  Stutzman 
and  Dishman  (1982),  can  be  assumed  to  fill  the  entire  path  homogeneously.  Accord¬ 
ing  to  that  model,  the  path-averaged  rain  rate  is 

R  =  R  when  R  <_  10,  and  is  for  R  >  10  given  by 

(30) 

'  '  R  ' 

For  slant  path  calculations,  the  rain  cell  extent  is  limited  to  height  level  h^ 

(0°  C)  and  the  secant  law  (13)  is  applied  for  predicting  elevation  angle  (i|») 
dependences. 
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Figure  3.  Specific  rain  attenuation  a  over  a  frequency  range  from 

1  to  1000  GHz  for  various  rainfall  rates,  R  =  0.25  to  150  mm/h 
using  a  four-segment  approximation  for  equation  (29). 


'  rc.  y.  . » 


The  MPM  is  restricted  to  approximate  only  the  most  fundamental  effects  by  rain, 
formulated  by  (28)  to  (30).  Considerable  refinements  need  to  be  added  if  other 
precipitation-induced  signal  degradations  (e.g.,  depolarization,  scatter  interfer¬ 
ence  between  links  sharing  the  same  frequency,  etc.)  are  to  be  predicted 
(Oguchi,  1983). 

3.  ABSORPTION  IN  THE  UPPER  ATMOSPHERE  (30  TO  100  km) 

An  auxiliary  model  to  the  MPM  calculates  specific  attenuation  due  to  isolated 
lines  above  h  =  30  km.  All  line  overlap  has  disappeared  and  even  closely-spaced 
doublets  in  the  0£  spectrum  (Table  1)  can  be  treated  as  isolated  lines.  Increasing 
height  implies  exponentially  decreasing  pressure  p,  which  causes  a  single  line  to 
narrow.  The  attenuation  peak  at  line  center  begins  to  diminish  when  the  absorber 
abundance  disappears  or  when  Equation  (20)  approaches  the  finite  Doppler  width 

yD  =  6.20vo/^m6  kHz  (31) 

where  m  is  the  molecular  weight  (e.g.,  O2,  m  =  32,  0  =  1 ,  vQ  =  60  GHz,  then 
Yd  =  65.8  kHz).  An  adequate  approximation  can  be  made  by  replacing  y  with 


while  retaining  the  first  term  (Lorentzian)  of  the  line  profiles  F"  (17)  and  F1  (18). 
Actually,  during  this  process,  the  line  profile  also  changes  from  a  Lorentzian  to  a 
Gaussian.  This  transition  is  described  by  a  Voigt  profile  (Liebe,  1981)  that,  for 
the  purpose  of  the  MPM,  can  be  passed  over.  Water  vapor  height  profiles  indicate 
that  e(h  >  20  km)  -  0,  except  for  a  very  small  residual .which  can  exist  up  to  about 
70  km.  Mesospheric  water  vapor  lines  can  be  neglected  for  the  MPM.  The  oxygen 
mixing  ratio,  however,  remains  constant  up  to  about  80  km  where  dissociation  to 
atomic  oxygen  begins.  At  120  km,  the  ratio  O2/O  has  dropped  off  to  0.25.  Oxygen 
lines  have  to  be  traced  above  30  km  and  the  situation  grows  complicated  due  to 
Zeeman  splitting  caused  by  the  geomagnetic  field. 

3.1  Oxygen  Zeeman  Patterns 

Zeeman  splitting  of  isolated  oxygen  lines  due  to  the  influence  ov  the  earth's 
magnetic  field  strength  H  (in  Gauss)  introduces  considerable  complications.  Each 
line  is  identified  by  a  quantum  number  K.  The  presence  of  the  steady  field  H  splits 
K1  lines  into  three  groups  of  (2^  +  1)  sublines,  thereby  redistributing  line  atten- 


uation  over  a  fixed  frequency  range.  In  principle,  three  different  Zeeman  patterns, 
al  ,2,3’  are  possible  for  any  K~Jine.  These  patterns  are  obtained  by  evaluating 
n  and  £  coefficients  in  tt  and  a1  groups.  Each  sum  of  sublines  determines  a  Zeeman 
pattern  according  to 

al  ,2,3  =  “o^/O  +  W  -  v|)/yh]2)  (33) 

where  the  scale  factor  aQ  =  a(vQ)  follows  with  (2),  (8),  (19),  (17),  (20),  (32) 
to  be 


ao  =  0J82  Vo(a1p/Yh)03exp[a2(1  '  0)^  dB/km  •  (34) 

The  line  center  frequency  of  each  Zeeman  component  is  given  by 

Vq  =  vQ  +  n  2.803  H  x  10'3  GHz  .  (35j 

A  scheme  to  calculate  number,  position  (n) »  and  intensity  (£)  of  individual  Zeeman 
components  was  derived  from  the  work  of  Lenoir  (1968)  by  Liebe  (1981;  Table  6).  The 
component  number  increases  with  quantum  number  r  (line  ID,  Table  1).  Numerical 
and  graphical  results  for  the  first  eight  lines  K1  =  1  to  7  illustrate  in  Table  5 
and  Figure  4  the  organization  of  Zeeman  components  into  three  groups  labeled 
tt,  a  ,  and  a  .  Each  component  is  treated  as  a  separate  independent  line  having 
identical  shape  and  width.  The  intensity  sum,  both  for  the  w  components  and  for  the 
o+o  components,  is  equal  to  the  integrated  intensity  of  the  unsplit  line  repre¬ 
sented  by  the  scale  factor  aQ  (34). 

The  three  main  patterns  a-j  (tt) ,a2(a+)  and  a 3(0")  for  lines  K*  =  1  to  29  (Table  1) 
have  been  calculated  at  local  p  -  0  -  H  conditions  of  the  U.  S.  Standard  Atmosphere 
(NOAA,  1976)  over  the  height  range  30  to  100  km  at  two  magnetic  field  strengths, 

H  =  0.3  and  0.6  G.  A  catalog  of  the  complete  set  is  presented  in  graphical  form  in 
Appendix  A.  Above  h  =  70  km,  individual,  now  mostly  Doppler-broadened  components 
become  discernible.  The  results,  in  combination  with  beginning  dissociation  to 
atomic  oxygen,  establish  h  =  100  km  as  the  boundary  to  outer  space  for  radio  path 
modeling.  Over  the  height  range  30  to  100  km  and  a  narrow  frequency  band  around 
vo  one  is  confronted  with  direction-dependent  (unisotropic)  attenuation  rates.  The 
first  step  towards  tying  the  three  main  Zeeman  patterns  to  a  global  coordinate 
system  for  H  and  antenna  position  and  expressing  them  for  two  principal  directions 
and  various  wave  polarizations  is  discussed  in  Section  4.2. 


*  *  ^  - 


24 


Table  5.  The  Frequency  Shifts  »7(K,M)  and  Relative  Intensity  Factors  f  (K,M) 
_ for  the  ZEEMAN  Components  of  the  K=l±  to  7±  Oxygen  Lines  (Table  i). 


^v.v.NV'V'VV.V.VW^.v.nv, 


.<  .v>- *.  v’srv*.^.  ■ \h^.T\Hvvv7' 


3.2  Trace  Gas  Spectra 

Atmospheric  trace  gases  such  as  carbon  monoxide  (CO),  ozone  (O^),  nitrous 
oxide  (N2O),  nitric  oxide  (NO),  and  sulfur  dioxide  (S02)  are  potential  millimeter 
wave  absorbers.  In  particular,  the  asymmetric  rotor  molecules  0^  and  SG?  exhibit  a 
rich  millimeter  wave  line  spectrum  (Rothman  et  al.,  1983b;  Park  et  al.,  1981). 
Scanty  abundances  (<  2  ppm/vol.)  of  trace  gas  absorbers  and  dominant  pressure¬ 
broadening  by  the  major  constituents  (=  p)  make  none  of  these  lines  detectable 
in  the  troposphere.  When  linewidths  (20)  have  narrowed  to  a  few  megahertz 
(h  _>  20  km),  stronger  lines  may  affect  the  millimeter  wave  behavior  of  radio  paths 
to  outer  space,  especially  if  abundances  have  increased  with  height,  as  is  the  case 
for  ozone.  Tiny  absorption  spikes  may  be  observed  on  an  earth-space  path  operating 
at  low  elevation  angles  in  transparent  portions  of  the  spectrum. 

Spectra  are  modeled  with  the  formulations  described  in  Section  2.1.1.  Tables 
of  the  spectroscopic  coefficients  have  to  be  prepared  before  line  spectra  of  trace 
gas  molecules  can  be  added  to  (8).  Several  millimeter  wave  lines  of  the  species 
CO  and  0^  have  been  observed  in  the  atmosphere.  Their  calculation  scheme  is 
detailed  below; 


CO  lines  in  air 

03  lines  in  air 

Equation 

ID 

c 

d 

(7), (8) 

vq,GHz 

Table  6 

Table  6 

S.kHz 

c1M(,10'1882exp[c2(l  -  0)] 

d1Md10“18a2*5exp[d2(l  -  0)] 

(36) 

Y.GHz 

c3p0°*8 

d3p0°'5 

(37) 

6 

0 

0 

Line  strengths  S  are  defined  to  be  proportional  to  the  number  density  M  in 
moleeules/cm  .  The  spectroscopic  coefficients  c^  to  Cy  d^  to  d^  are  given  in 
Table  6.  Models  exist  for  the  height  profiles  M  H(h)  (N0AA,  1976).  For  example, 
ozone  displays  a  maximum  •-  5  *  10  Ccm"j  in  the  profile  N^(h)  at  a  height  of 
h  =  22  km  (p  =  4  kPa,  0  =  1.37)  where  the  air  density  is  Na  -  1.3  *  1018cm"3. 

Under  these  conditions,  the  peak  attenuation  is  u(vQ)  3  7.6  *  10“**  dB/km  for  the  line 
centered  at  =  101.7  GHz  (Table  6). 


Table  6.  Carbon  monoxide  (CO)  and  selected  Ozone  (0^) 
line  parameters  (Waters,  1976;  Rothman  et  al.,  1981) 
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In  principle,  line  tables  for  trace  gases  could  be  included  in  the  MPM;  in 
practice,  it  turns  out  that  these  lines  are  relatively  weak,  that  their  main 
importance  lies  more  in  remote  sensing  applications  of  abundances  and,  most  impor¬ 
tant,  that  their  assessment  requires  considerable  computing  time.  For  these 
reasons,  so  far,  we  have  omitted  trace  gas  spectra  from  the  MPM. 

4.  MODELING  STRATEGY 

The  MPM  systematizes  radio  path  parameter  calculations  (see  Section  1.1)  in 
two  computer  routines  PI  and  P2.  Routine  PI  constitutes  the  main  model,  while  the 
auxiliary  routine  P2  is  activated  only  if  line  details  over  the  height  range 
h  =  30  to  100  km  are  required.  This  simplifies  the  handling  of  the  large  number 
of  spectroscopic  quantities  and  computational  steps  involved.  The  computations 
are  performed  in  a  user-friendly,  interactive  mode  designed  to  minimize  possibilities 
for  error.  An  operator  needs  to  specify  a  frequency  range  f,(>l  GHz)  to 
^2 (— ^ 0^0  GHz)  and  the  following. 

Horizontal  path:  height  level  h,  model  atmosphere  or  set  of  appropriate 
p,T-values,  relative  humidity  RH,  hydrosol  concentration  w, 
rainfall  rate  R,  and  path  length  L. 

Slant  path:  initial  height  level  hQ,  elevation  angle  ip,  relative 

humidity  RH  at  hQ,  final  height  h^.,  and  model  atmosphere. 

For  an  earth-space  path  h^  °  30  km,  unless  routine  P2  is 
called  to  include  oxygen  line  details  (vQ  +  50  MHz)  up  to 
hf  =  100  km. 

Routine  PI  generates  a  frequency  grid  with  a  resolution  tied  to  the  line  center  fre¬ 
quencies  that  fall  within  the  desired  range.  Two  consecutive  values,  vQ  and  v^, 
are  taken  and  the  number  of  points  between  them  is  specified  (usually  2  to  5).  Also 
automatically  added  to  the  grid  are  the  half-power  points  at  vQ  ♦  y.  Thus,  it  is 
assured  that  attenuation  maxima  and  dispersion  peaks  will  be  plotted.  The  program 
feeds  one  frequency  at  a  time  into  the  subroutine  and  receives  values  of  u  and  tS, 
which  are  stored  temporarily  in  two  files.  When  all  numbers  of  the  grid  have  been 
handled,  a  numerical  integration  is  performed  layer-by-layer  over  the  height  range 
of  the  desired  radio  path  and  the  output  is  presented  in  a  numerical  printout  and 
processed  by  a  graphics  routine.  In  the  troposphere  (h  <  10  km),  only  absorption 
by  HgO  (vapor  and  liquid)  and  Oj  is  important.  Above  h  *  5  km,  spectral  lines  of 


weak  trace  gas  abundances  (0^,  CO,  ^0,  etc.)  become  detectable  over  narrow 
(v  +  0.2  GHz)  frequency  ranges  (Waters,  1976;  Poynter  and  Pickett,  1981),  but  they 
are  not  included  in  the  MPM,  in  keeping  with  the  goal  of  modeling  simplicity.  Also, 
computational  economy  cannot  be  overlooked,  since  each  specific  value  of  a  and  8 
requires  the  summation  of  up  to  82  contributions  (8)  at  one  frequency  and  the 
height  integration  typically  is  carried  out  over  as  many  as  56  layers. 

4.1  Spectroscopic  Data  File  and  Synthetic  Atmospheres 
The  system  of  attenuation/dispersion  routines  consists  of  four  FORTRAN  4  pro¬ 
grams  and  11  permanent  data  files.  Of  these  data  files.  No.  1  is  the  permanent 
spectroscopic  line  data  file  LNDAT  which  contains  line  center  frequencies  and  other 
spectroscopic  parameters  (strength,  width,  interference)  for  each  of  48  0^  lines 
and  30  H20  lines  in  air  (Tables  1  and  2).  Files  No.  2  to  11  are  synthetic  atmos¬ 
phere  data  bases  that  list  height,  pressure,  temperature,  and  maximum  relative 
humidity  data  for  45  discrete  height  levels  (Valley,  1965);  they  are  called  by  the 
following  names: 

STDATM  1  -  U.  S.  Standard  Atmosphere  1976, 

STDATM  2  -  Tropical  15  degrees  north, 

STDATM  3  -  Subtropical  30  degrees  north,  July, 

STDATM  4  -  Subtropical  30  degrees  north,  January, 

STDATM  5  -  Midlatitude  45  degrees  north,  July, 

STDATM  5  -  Midlatitude  45  degrees  north,  January, 

STDATM  7  -  Subarctic  60  degrees  north,  July, 

STDATM  8  -  Subarctic  60  degrees  north,  January, 

STDATM  9  -  Arctic  75  degrees  north,  July, 

STDATM  0  -  Arctic  75  degrees  north,  January. 

The  permanent  line  data  file  and  a  model  atmosphere  are  assigned  to  each  job  as 
two  local  file  tapes  (TAPE6,7). 

4.2  PI:  Lower  Atmosphere  Program  (0  to  30  km) 

There  are  three  FORTRAN  4  programs  that  make  up  the  attenuation/delay  system 
for  this  height  range.  The  file  names  are  as  follows: 
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SETPLT  -  interactive  run  set-up  program; 

ATNPLT  -  computation  program  according  to  the  physical  model 

for  N,  including  the  revisions  discussed  in  Section  2; 

ATNOUT  -  output  routine  producing  listings,  plots,  and 

numerical  integration  for  cumulative  calculations 
through  the  array  of  discrete  height  levels. 

A  control  language  procedure  file  DERP  is  used  to  properly  set  up  a  job. 

The  three  programs  communicate  with  each  other  through  several  local  files. 

Local  files  TAPE 1  to  TAPE10  are  reserved  for  these  programs.  The  first  program  is 
interactive  and  requests  from  the  user  the  frequency  range,  climatological  data, 
plot  sizes,  and  several  output  options.  It  then  writes  the  input  information  to 
local  file  TAPE1,  the  line  data  base  used  for  calculations  is  assigned  to  TAPE2, 
and  if  a  plot  is  requested,  the  plot  parameters  are  stored  in  TAPE5. 

The  second  program,  in  generating  the  attenuation/delay  data,  performs  the 
majority  of  all  calculations.  This  program  sets  up  a  grid  of  frequency  (x-axis) 
values  consisting  of  any  combination  of  line  center  frequencies,  half-power  points, 
and  even  increment  values.  Then,  if  a  plot  is  desired,  the  program  automatically 
adds  additional  resolution  points  to  the  frequency  grid  so  that  a  smooth  curve  is 
obtained.  On  cumulative  path  calculations,  no  half-power  and  additional  resolution 
points  are  added  because  this  mode  requires  considerably  more  computational  effort. 

In  this  case,  the  grid  increment  has  to  be  chosen  for  the  desired  frequency  range 
as  a  compromise  between  cost  and  needed  resolution.  The  elements  of  the  frequency 
grid  are  then  passed  step-by-step  through  the  propagation  model,  and  two  arrays  with 
attenuation-vs. -frequency  and  delay-vs. -frequency  data  for  each  set  of  climatological 
data  are  completed,  each  constituting  a  “curve."  The  points  are  written  to  binary 
local  files  TAPE3  and  TAPE4,  curve  by  curve. 

The  third  program  handles  output.  Based  on  user  specifications,  a  listing  and 
a  pM  (if  desired)  are  produced.  The  listing  is  sent  to  local  file  1APE1G,  and  the 
plots  are  made  using  the  DISSPIA  plotting  routine  In  conjunction  with  the  FR-80 
microfilm  plotter.  If  a  horizontal  path  mode  was  specified,  the  output  program 
prints  a  heading  that  consists  of  climatological  data  and  hQ.  Multiple  curves  are 
combined  in  the  same  plot  frame.  If  a  cumulative  (zenith)  path  was  specified,  the 
output  prints  a  table  of  climatological  data  used  at  each  height  level.  Attenuation 
and  delay  data  for  these  levels  are  summed  to  approximate  integration  over  height. 


For  each  frequency,  attenuation  and  delay  values  on  two  consecutive  height  levels 
are  averaged,  then  multiplied  by  the  height  difference  (layer  thickness  within 
model  atmosphere).  Zenith  attenuation  (dB)  and  delay  (s>  are  accumulated  at 
each  frequency  as  the  program  moves  from  the  specified  initial  to  the  final  height 
level.  Large  frequency  ranges  and  small  frequency  grid  elements  should  be  avoided 
to  keep  expenses  under  control. 

Details  for  user-specified  parameters  are  discussed  below. 

1.  RUN  TYPE  —  a  code  number  1  to  6  specifies  the  following  outputs: 

EVEN  GRID  LINE  CENTERS 
no  yes 

yes  no 

yes  yes 

no  yes 

yes  no 

yes  yes 

A  repeat  run  may  be  done  by  specifying  a  negative  run  type  code.  For  example, 
run  code  3  is  specified,  then  the  user  decides  he  would  like  a  plot  of  the  same 
data,  By  rerunning  the  system  and  specifying  a  code  (-4),  a  plot  is  produced  using 
identical  climatological  data  but  without  even  grid  frequencies.  This  is  cheaper 
and  produces  a  neater  output  than  specifying  code  6,  since  the  listing  does  not  eon 
tain  the  resolution  points  used  in  plotting. 

2.  FREQUENCY  RANGE  --  beginning  and  ending  frequency  are  specified.  For 
calculations  from  30  to  300  GHz,  the  proper  response  is  “30.,  300.“. 

3.  EVEN  GRID  --  this  number  is  the  increment  for  an  even  frequency  grid.  For 
2  GHz  steps,  the  response  is:  “2.“ 

4.  WATER  VAPOR  CONTINUUM  ABSORPTION  -  is  modeled  by  the  expression 

a  (bfepfV  ♦  bse2fxV‘)10*6  ppm 

Coefficients  bf  $  and  exponents  x,y  can  be  specified,  or  the  default  values  (new 
continuum  model ,  Eq.  24) 


K. 


EVu 


bf  =  1.40  ,  x  =  1  ,  y  =  2.5 

b$  =  54.1  ,  x'  =  1  ,  y*  =  3.5 

are  used  automatically. 

5.  LINE  STRENGTH  CUT-OFF  —  in  order  to  reduce  computer  run  time,  a  cut-off  value 
can  be  applied  to  the  resonance  line  strength  factors.  All  lines  in  the  permanent 
file  are  checked,  and  those  falling  below  the  specified  cut-off  are  disregarded  in 
attenuation/delay  calculations. 

6.  PATH  MODE  --  horizontal  or  zenith  paths  may  be  selected.  In  the  horizontal 
path  mode,  discrete  height  h,  pressure  p,  temperature  T,  relative  humidity  RH,  and 
droplet  concentration  w  or  rain  rate  R  for  each  curve  are  entered.  Combinations 
of  h(p,T)  can  be  also  drawn  from  the  10  atmospheric  local  files. 


ft 


In  the  slant  path  mode,  a  model  atmosphere  and  elevation  angle  ^  are  specified, 


together  with  a  height  range  between  0  and  30  km  and  a  relative  humidity  at  hQ.  A 


cloud  height  range  may  be  defined  when  100  percent  relative  humidity  conditions 
exist,  where  the  mean  droplet  concentration  w  is  entered.  For  example,  "2.,  4.,  .1“ 
establishes  a  cloud  between  2  and  4  km  with  a  droplet  concentration  of  w  =  0.1  g/m3. 
Zenith  path  behavior  is  calculated  when  ^  =  90°. 


7.  PLOT  PARAMETERS  --  if  a  plot  has  been  requested,  the  program  asks  for  plot 
axis  length  and  labeling  inputs. 


8.  ERROR  CORRECTION  —  at  this  point,  if  any  errors  in  entry  have  been  made,  or 
if  the  user  changes  his  mind,  any  input  parameter  may  be  changed.  A  final  reply 
"NO"  to  the  option  of  initiating  changes  executes  the  program  and  produces  the  re¬ 
quested  output  in  microfilm  plots  and/or  printed  listings  on  local  file  TAPE10. 


4.3  P2:  Isolated  Lines  in  the  Upper  Atmosphere  OQ  to  100  km) 

The  routine  P2  calculates  the  response  of  isolated  lines  over  the  height  range 
h  =  30  to  100  km  as  formulated  in  Section  3.1.  The  routine  takes  into  account  the 
complicated  Zeeman  effect  of  the  0£  microwave  lines.  Formally,  each  oxygen  splits 
proportionally  with  the  geomagnetic  field  strength  H  Into  numerous  sub-lines,  which 
superimpose  to  three  main  Zeeman  patterns  ^  ^  spread  over  a  megahertz  frequency 


33 


scale  (see  Appendix  A).  Routine  P2  stops  at  this  point  and  additional  information 
has  to  be  provided  before  a  mesospheric  02  line  pattern  az  can  be  calculated. 

The  propagation  medium  is  anisotropic  in  close  vicinity  of  02  spectral  lines. 
Transformation  formulae  which  reduce  the  three  Zeeman  components  2  3  into  an 
observable  attenuation  rate  az  with  respect  to  given  coordinates  for  the  geomagnetic 
field  vector  and  the  propagating  radio  wave  (or  receiving  antenna)  have  been  given 
by  Hartmann  and  Kunzi  (1983).  Here,  only  the  simpler  case  of  angle  dependence  be¬ 
tween  wave  plane  and  geomagnetic  field  direction  is  discussed.  For  linearly 
polarized  radiation,  an  angle  of  orientation  <J>  is  defined  between  the  wave's 
magnetic  field  component  in  a  plane  of  constant  phase  (i.e.,  perpendicular  to  the 
direction  of  propagation)  and  the  geomagnetic  field  direction  H.  For  circularly 
polarized  radiation,  the  angle  <f>  is  defined  between  the  plane  of  constant  phase  and 
H.  A  mesospheric  02  line  az  then  generally  consists  of  a  mixture  of  the  three 
Zeeman  patterns: 

Wave  polarization  Zeeman  line  a7  = 

2  2 

linear  0^  cos  <J>  +  (a2  +  a^)  sin  <j> 

2  2 

circular,  right-handed  2a2sin  <J>  +  0.5(a^  +  a2  +  a3)cos  <J>  (38) 

2  2 

circular,  left-handed  2a3sin  $  +  0.5(a1  +  a2  +  a3)cos  <p  . 

These  equations  assure  that  in  all  cases  the  sum  of  the  Zeeman  components  equals  the 
integrated  absorption  strength  of  the  unsplit  (H  =  0)  line. 

5.  RESULTS 

The  quest  for  a  model  that  accurately  describes  the  mean  of  millimeter  wave 
properties  of  the  atmosphere  in  terms  of  measurable  primary  parameters  has  been  a 
goal  of  the  ITS  millimeter  wave  program.  The  current  status  of  an  MPM  (ITS  version) 
has  been  described  in  the  preceding  sections.  The  model  provides  a  capability  to 
forecast  atmospheric  transmission  effects  on  communications  and  sensor  systems. 

The  influence  of  the  atmosphere  was  structured  in  four  categories;  namely,  dry  air 
(p,u),  humidity  (v.RH),  hydrosol  (w),  and  rain  (R)  for  the  lower  part,  h  =  0  to  30 
km,  and  dry  air  (p,6 )  and  magnetic  field  strength  (H)  for  the  upper  part,  h  =  30  to 
100  km.  In  this  section,  graphical  and  numerical  results  are  presented  in  the 
form  of  predictions  over  a  frequency  range  from  10  to  300  GHz  by  assuming  typical 


climatic  conditions.  Part  of  the  results  can  be  validated  by  comparison  with  data 
from  field  and  laboratory  measurements,  of  which  there  are  presently  but  a  few  above 
100  GHz;  another  part  (i.e.,  the  opaque  medium,  A  >  40  dB),  will  always  rest  on  the 
underlying  theory.  As  more  data  become  available,  further  revisions  and  extensions 
of  the  MPM's  data  base  will  have  to  be  considered. 

5.1  Horizontal  Radio  Paths 

A  homogeneous  absorber  population  characterizes  a  horizontal  radio  path.  One 
set  of  p-0-v  +  w  +  R  combinations  is  used  to  predict  the  specific  rates  of  a  and  3. 
The  dispersive  delay  3(f)  is  presented  in  three  units  of  measure  (see  Equation  3): 
3(f)  in  radians/km,  3*(f)  in  ns/km,  or  D(f)  in  ppm.  A  typical  output  is  illustrated 
in  Figure  5.  Dry  and  moderately  humid  air  (curves  1,2)  under  sea  level  conditions 
allow  four  millimeter  wave  window  ranges  W1  to  W4.  Transmission  through  these 
windows  is  reduced  by  high  humidity,  fog,  and  light  rain  (curves  3  to  5),  and  is 
seriously  impeded  by  steady  rain  end  thundershowers  (curves  6,7).  The  delay 
spectrum  3*  from  gaseous  resonance  contributions  displays  frequency  gradients 
around  the  line  centers  that  can  distort  broadband  signal  transmission  (e.g., 
ns-pulses).  Details  around  the  22  GHz  water  vapor  line  are  enhanced  in  Figure  6. 

This  relatively  weak  line  imposes  no  serious  effects  upon  signal  transmissions.  As 
demonstrated  in  Figure  7,  the  influence  of  water  vapor  becomes  important  above 
70  GHz  and,  with  increasing  frequency,  actually  dominates  more  and  more  the  milli¬ 
meter  wave  properties  of  clear  air. 

The  frequency  range  45  to  75  GHz  is  dominated  by  oxygen  absorption,  as  evident 
from  Figure  8.  The  change  in  attenuation  from  h  =  0  to  6  km  is  small.  The  lines 
are  merged  into  a  more  or  less  unstructured  band  with  a  pressure-proportional  atten¬ 
uation  maximum.  Above  h  «  10  km,  however,  one  can  see  a  significant  difference--the 
unstructured  60  GHz  band  shape  transforms  into  a  highly-structured  isolated  line 
response  with  distinct  peaks  and  valloys.  The  lines  narrow  with  increasing  altitude, 
deepening  the  valleys  while  the  line  peaks  stay  approximately  constant.  This 
behavior  is  detailed  in  Appendix  A. 

Tracing  the  line  peaks  to  zero  intensity  over  the  range  h  =  30  to  100  km  intro¬ 
duces  complications  of  Zeeman  splitting  by  the  geomagnetic  field  H,  which  has  been 
treated  with  program  P2.  Model  calculations  for  0,,  lines  labeled  1C  =  1  to  29  (see 
Table  1)  are  displayed  in  Appendix  A  assuming  two  geomagnetic  field  strengths, 

H  -  0.3  and  0.6  Gauss.  Numerical  values  for  peak  attenuations  of  the  a-|  (7r)-patterns 
are  listed  in  Table  7.  A  comparison  between  programs  PI  and  P2  at  h  =  30  km  for  the 
unsplit  line  (H  =  0)  peak  reveals  the  validity  of  switching  at  that  point  to  a 
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Figure  6.  Specific  attenuation  a  and  refractive  dispersion  D  over  a 
frequency  range  from  0.5  to  50  GHz  for  saturated  (RH  =  100%) 
sea  level  air  (p  =  101.3  kPa)  at  four  temperatures. 


Table  7.  Peak  attenuation  calculated  for  the  Zeeman  it  patterns  of 
K~  =  1  to  29  with  Program  P2  for  magnetic  field  strengths 
H  =  0,  0.3,  and  0.6  6  assuming  a  U.S.  Standard  Atmosphere 
(NOAA,  1976). 

a)  Specific  attenuation  a-j(v  )  for  h  =  30  (see  also  PI  for 

comparison),  40,  50,  _ °100  km  (see  Appendix  A). 

b)  Zenith  attenuation  A^(vq)  for  h  =  30  to  100  km. 
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the  unsplit  line  (H  =  0)  peak  reveals  the  validity  of  switching  at  that  point  to  a 
single-line  treatment.  The  error  is  less  than  -6  percent  for  the  four  doublets  and 
less  than  -1  percent  for  lines. 

Line  intensities  are  sensitive  to  temperature  variations  at  various  rates. 
Useful  approximations  at  a  fixed  frequency  are  the  power  laws 

a(T)  =  a(3OOK)0y  and  0(T)  =  6(3OOK)0y'  .  (39) 

For  sea  level  air,  we  obtain 


f.GHz 

50 

55 

60 

65 

70 

119 

120 

y 

2.49 

1.30 

2.35 

1.20 

2.72 

2.13 

2.70 

y* 

1.84 

1.98 

3.30 

1.87 

1.88 

1.96 

1.98 

5.2  Zenith  Path  Behavior 

Cumulative  transfer  characteristics  of  the  neutral  atmosphere  are  evaluated  by 
a  numerical  integration  of  (9)  and  (10).  For  a  zenith  path,  the  increments  dx  are 
equal  to  dh,  and  dh  is  chosen  so  that  each  slab  is  quasi -homogeneous  (p,T,v,w  and  R 
are  constant).  The  formalism  developed  in  Sections  2  to  4  is  applied  to  calculate 
a(dh)  and  0(dh).  Height  profiles  of  the  meteorological  quantities  are  obtained  from 
model  atmospheres.  A  standard  example  is  the  one-way  zenith  response  through  the 
U.  S.  Standard  Atmosphere  (NOAA,  1976).  The  calculation  is  performed  with  PI 
u layer- by- layer, “  encompassing  48  slabs  between  h  =  0  and  30  km.  Examples  for 
zenith  attenuation  and  delay  are  presented  in  Figures  9  to  11.  One  notices  that 
outside  the  0^  absorption  frequency  ranges  (50  to  70  GHz  and  115  to  123  GHz),  water 
vapor  content  V  (15)  and  cloud  water  content  W  (16)  play  a  dominant  role.  The 
refractive  delay  of  the  air  mass  (14)  is  B*  °  7.62  ♦  0.09  *  0.21  V  ns,  where  the 
0.09  ns  term  is  the  contribution  from  the  height  range  30  to  100  km.  Selected  atten¬ 
uation  values  Az  are  collected  in  Table  8. 

One  can  attempt  to  prediet  zenith  path  attenuation  A^  from  the  specific  atten¬ 
uation  rate  aQ(h0)  by  reducing  the  inhomogeneous  atmosphere  to  an  equivalent  homo¬ 
geneous  “thiekness“ 
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Figure  10.  One-way  zenith  phase  dispersion  Bz  over  a  frequency  range  from 
1  to  300  GHz  for  three  values  of  water  vapor  content  V  (15)  ov*r 
the  height  range  h  =  0  to  30  km  (PI)  assuming  a  cloud-free  0.  $. 
Standard  Atmosphere. 
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Figure  11. 


One-way  zenith  path  delay  B*  and  attenuation  A  over  a  frequency 
range  from  15  to  25  GHz  tiwugh  the  eloud-free2U.  S.  Std.  Atm.  76 
for  live  values  of  water  vapor  content  V  over  the  height  range 
h  =  0  to  30  km  (PI). 
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Table  8.  One-way  zenith  attenuation  A  through  the  U.  S.  Standard  Atmosphere 
at  frequencies  dominated  by  water  vapor  absorption,  covering  the 
range  20  to  300  GHz. 


MODEL  PI 

Frequency  Az(50%)*  Az(100%)* 

a(30) 

f  (h  =  0-30  '■  n) 

(h  =  30) 

GHz 

08 

dB/kra 

20 

.27 

.44 

.00 

22. 235(^0) 

.54 

.34 

.91 

.57 

.00 

.00 

30 

.24 

.35 

.00 

32 

.25 

.25 

.00 

35 

.28 

.39 

.00 

40 

.40 

.52 

.00 

45 

.68 

.82 

.00 

50 

1.64 

1.81 

.00 

60 

Oxygen  Line  Absorption 

70 

2.27 

2.66 

.00 

80 

.96 

1.37 

.00 

90 

.37 

1.38 

.00 

100 

.97 

1.60 

.00 

105 

1.08 

1.78 

.00 

no 

1.33 

2.10 

.00 

115 

2.57 

3.42 

.00 

120 

9.22 

10.1 

.00 

125 

1.78 

2.79 

.00 

130 

1.53 

2.64 

.00 

no 

2.45 

2.99 

.00 

150 

1.99 

3.62 

.00 

160 

2.59 

4.75 

.00 

170 

4.46 

8.25 

.00 

160 

30.2 

55.6 

.00 

183.3»{H.O) 

135  « 

93.3 

60.7 

168 

109 

.11 

.00 

190 

13.1 

24.3 

.00 

200 

5.00 

9.26 

.00 

210 

4.14 

7.63 

.00 

220 

4.06 

7.52 

.00 

230 

4.25 

7.82 

.00 

240 

4.52 

8.12 

.00 

250 

4.86 

8.93 

.00 

260 

5.24 

9.64 

.00 

270 

5.17 

10.4 

.00 

280 

6.13 

11.4 

.00 

290 

6.82 

12.6 

.00 

300 

7.74 

14.3 

.00 

•The  cumulative  values  A  { 50t ) ,  A  (lOtti)  are  for  water  vapor  contents  V  ■  14.4  ai 
28.7  m,  respectively;  *(30)  is  the  specific  attenuation  at  the  level  h  «  30  km 
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The  relations  between  A  and  a  have  been  calculated  with  PI  and  examples  for 

z  o 

hQ  =  0  km  and  the  U.  S.  Standard  Atmosphere  are  listed  below  (see  also  Figure  11): 


f,GHz 

15 

22 

35 

60* 

90 

119* 

140 

183 

220 

V,mm 

h$,km 

0 

5.0 

5.0 

5.0 

15.4 

4.8 

50 

4.8 

4.5 

4.6 

15 

3.1 

3.0 

2.9 

15.3 

2.2 

35 

1.9 

3.2 

1.8 

30 

2.4 

2.9 

2.3 

15.2 

1.9 

25 

1.7 

3.2 

1.7 

60 

2.2 

2.6 

2.1 

15.1 

1.8 

20 

1.7 

3.2 

1.7 

★ 

Valid  for  h  =  0  to  30  km;  see  Table  9  for  contributions 
over  the  range  h  =  30  to  100  km. 


Scale  heights  hg  assume  different  values  depending  on  frequency  f  and  water  vapor 
content  V,  which  makes  the  effective  height  concept  not  very  universal. 

Phase-coherent  satellite  beacon  signals  at  frequencies  around  the  22  GHz  line 
experience,  in  principle,  different  delays  Bz  due  to  a  line  component  which  depends 
linearly  upon  water  vapor  content  V.  This  is  demonstrated  in  Figure  11.  Analysis 
reveals  that  the  water  vapor  group  delay  between  two  signals  placed  symmetrically 
about  vQ  is  extremely  small;  for  example,  when 

a)  frf2  =  \>Q  +  0.5  GHz:  Bf  -  *  (0.0087  V  -  0.07)  ps  and 

Az  « (0.07  +  C  028  V)  dB;  or 

b)  frf2  =  \>Q  +  2.0  GHz:  Bf  -  B£  *(0.020  V  -  0.27)(max.)  ps  and 

Az  *{0.07  +  0.017  V)  dB. 

The  frequency-dependent  delay  (14)  of  each  signal  traversing  the  moist  air  mass  is 
B*  ■  (7747  +  22.1  V)  ps.  The  22  GHz  water  vapor  line  is  not  a  good  water  vapor 
detector.  Sensitivity  can  be  increased  by  a  factor  of  3  by  looking  at  the  signals 
under  an  elevation  angle  20°,  or  by  a  factor  of  about  180  by  shifting  to  the 
183  GHz  line. 

Examples  of  zenith  attenuation  through  a  cloudless  atmosphere  are  shown  for  the 
frequency  range  10  to  350  GHz  in  Figure  12.  The  atmosphere  is  opaque  for  any  system 
attempting  to  operate  from  the  ground  to  outer  space  when  Az  >  30  dB.  The  55  to 
65  GHz  band  has  especially  stable  shielding  qualities.  The  shielding  factor  of  the 
air  mass  has  to  be  evaluated  by  analytical  means;  the  correct  attenuation  at  the 
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oxygen  line  centers  calls  for  a  detailed  analysis  with  P2.  An  example  of  such  an 
exercise  is  exnibited  in  Figure  13.  The  numerical  account  of  the  integration  over 
h  =  30  to  100  km  at  selected  frequencies  can  be  found  in  Table  9  for  the  opaque 
region  (a)  and  for  semi-transparent  regions  (b)  to  either  side  of  the  60  GHz  band. 
Maximum  shielding  (A  150  dB)  of  satellite-based  systems  against  interference  from 
ground  levels  can  be  achieved  between  58.9  and  61.2  GHz. 

For  earth- to-space  radio  paths,  both  climatic  variations  in  the  state  of  the 
atmosphere  and  elevation  angle  dependence  determine  communications  link  performance 
at  a  given  frequency.  Assessing  the  low  elevation  behavior  of  A  by  evaluating  ray 
curvature  with  (12)  yields  the  following  selected  results: 

Cumulative  Path  Attenuation  A 


1 p  =  90° 

(zenith) 

=  10° 

(slanted)  ip 

=  0°  (tangential) 

f 

(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

GHz 

dB 

dB 

dB 

30 

0.13 

0.35 

0.75 

2.0 

7.4 

30 

70 

2.0 

2.7 

11.5 

15.3 

115 

175 

90 

0.31 

1.4 

1,8 

8.0 

19 

126 

140 

0.26 

3.0 

1.5 

17.2 

18 

270 

(a)  V  = 

1 .5  mm 

,  (b)  V  = 

30  mm 

(see  Figure  12, 

Table  8) 

This  section  has  presented  results  of  calculations  with  the  ITS  version  of  an 
MPM.  Atmospheric  transfer  properties  in  the  millimeter  wave  range  are  different 
from  those  in  the  lower  frequency  bands-- the  lower  air  mass  can  become  opaque  depend¬ 
ing  on  frequency,  path  length,  elevation  angle,  and  climatic  conditions. 

6.  SUMMARY  AND  RECOMMENDATIONS 

The  neutral  atmosphere  {h  <  100  km)  has  been  characterized  as  a  propagation 
medium  for  millimeter  waves.  The  basic  physics  of  signal  degradations  by  the  atmos¬ 
phere  has  been  cast  into  a  model  (MPM)  capable  of  predicting  attenuation  and  delay 
effects  caused  by  moist  air  containing  suspended  and  precipitating  water  drops.  The 
MPM  was  streamlined  for  optimum  computer  run  time  but  without  undue  approximations 
for  the  clear  >r  case.  Exact  line-by-line  calculations  are  performed  below  1000 
GHz  *or  78  resonance  features  of  the  gases  0£  and  H2O.  Experimentally  confirmed  dry 
and  net  continuum  absorption  terms  (22)-(24)  have  been  added  to  account  for  attenua¬ 
tion  in  between  lines.  Nevertheless,  more  than  300  spectroscopic  coefficients  are 
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Figure  13.  One-way  zenith  attenuation  A  over  a  frequency  range  from  53  to 
67  GHz  dominated  by  oxygen  Itne  absorption  for  V  =  15  mm  over  the 
height  range  h  *  0  to  100  km  (PI  +  P2).  The  line  peaks  are  for 
it  patterns  and  H  =  0.3  G  (see  Table  7). 
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Table  9.  One-way  zenith  attenuation  A^  through  the  U.S.  Standard  Atmosphere 
at  selected  frequencies  (max/min  values)  dominated  by  oxygen  absorp¬ 
tion:  (a)  Opaque  Region  (>  30  dB)  and  (b)  Semi-Transparent  Region 
(<  30  dB)  (see  Fig.  13). 

The  cumulative  values  A(50%),  A(1 00%)  are  for  V  =  14.4  and  28.7  mm, 
respectively;  a(30)  is  the  specific  attenuation  at  h  =  30  km;  A~  ~ 
and  Ag  g  are  contributions  of  mesopheric  Zeeman  fr-patterns  assuming 
a  geomagnetic  field  strength,  H  =  0.3  or  0.6  Gauss. 

(a)  Opaque  Region  (>  30  dB) 
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Digit  In  parenthesis  gives  the  Standard  Deviation  from  the  mean  (in  terms  of  the 
final  listed  digits)  due  to  magnetic  field  (H  =  0.3  to  0.6  G)  and  humidity 

(V  3  0  to  30  mm)  variations. 


(b)  Semi-Transparent  Region  (<  30  dB) 


Line 

Frequency  10 

f  K± 

MODEL  PI 

MODEL  P2 

MODEL  PI 

+  P2 

A( 50%) 
(h  =  0 

A(100X) 
-30  km) 

ct{30) 

(h  =  30) 

A0  3  (7r)  Aq  g 

(h  =  30  -100  km) 

Az 

(h  =  0  -100  km) 

GHz 

dB 

dB/km 

dB 

dB 

50.99 

35- 

2.19 

2.38 

.00 

.03 

.02 

2  3(1)* 

51.25 

2.40 

2.59 

.00 

2.5(1) 

51.50 

33- 

2.66 

2.86 

.00 

.07 

.06 

2.8(1) 

51.76 

2.96 

3.16 

.00 

3.KD 

52.02 

31- 

3.38 

3.58 

.00 

.18 

.15 

3.6(2) 

52.28 

3.82 

4.02 

.00 

3.9(1) 

52.54 

29- 

4.53 

4.74 

.01 

.46 

.38 

5.1(2) 

52.81 

5.20 

5.41 

.00 

5.3(1) 

53.07 

27- 

6.46 

6.68 

.02 

1.10 

.91 

7.6(2) 

53.33 

7.44 

7.66 

.00 

7.7(1) 

53.60 

25- 

9.75 

9.97 

.05 

2.45 

2.06 

11.9(3) 

53.86 

11.1 

11.3 

.00 

11.2(1) 

54.13 

23- 

15.3 

15.6 

.12 

5.09 

4.31 

20.2(4) 

54.40 

17.0 

17.2 

.00 

17.1(1) 

54.67 

21- 

24.5 

24.7 

.24 

9.85 

8.39 

33.7(8) 

54.95 

26.2 

26.4 

.00 

26.3(1) 

Oxygen  Band  Center  (Table  9a) 

64.95 

23.6 

24.0 

.00 

23.8(2) 

65.22 

23+ 

21.6 

21.9 

.17 

7.50 

6.40 

28.8(5) 

65.49 

15.7 

16.1 

.00 

15.9(2 

65.77 

25+ 

13.9 

14.3 

.08 

3.73 

3.18 

17.7(4) 

66.03 

10.7 

11.0 

.00 

10.8(2) 

66.30 

27+ 

9.35 

9.73 

.03 

1.73 

1.47 

11.11 

4) 

66.57 

7.53 

7.92 

.00 

7.71 

:2) 

66.84 

29+ 

6.62 

7.02 

.01 

.74 

.63 

7.5(3) 

67.10 

5.60 

6.00 

.00 

5.8| 

I2) 

67.37 

31+ 

4.99 

5.40 

.01 

.30 

.27 

5.51 

*2) 

67.64 

4.40 

4.81 

.00 

4.61 

(2) 

67.90 

33+ 

3.99 

4.40 

.00 

.12 

.10 

4.51 

[2) 

68.17 

3.62 

4.03 

.00 

3.8(2) 

68.43 

35+ 

3.33 

3.74 

.00 

.04 

.03 

3.5(2) 

Digit  in  parenthesis  gives  the  Standard  Deviation  from  the  mean  (in  terms  of  the 
final  listed  digits)  due  to  magnetic  field  (H  =  0.3  to  0.6  G)  and  humidity 
(V  5  0  to  30  mm)  variations. 


still  involved  to  various  degrees  in  the  calculation  of  one  output  value.  The  pre¬ 
sented  MPM  version  is  physically  sound,  has  been  tested  with  available  data  (see 
Appendix  C),  and  was  organized  to  adapt  easily  to  requirements  for  communications 
applications.  The  MPM,  which  inputs  routinely  measured  meteorological  parameters, 
can  serve  as  a  source  for  simpler,  engineering-type  approximations  that  suffice 
when  applications  are  limited  in  frequency  or  height  range  (Pierluissi  et  al.,  1982; 
Liebe,  1983). 

The  MPM  consists  of  a  main  program  PI  for  the  lower  atmosphere  (<  30  km)  sup¬ 
plemented  by  a  second  program  P2  for  the  upper  atmosphere  (<  100  km).  Both  are 
written  in  FORTRAN  IV  and  structured  to  be  operated  interactively  in  a  user-friendly 
mode  by  nonexperts.  Test  runs  provided  a  wide  selection  of  numerical  and  graphical 
examples,  which  elucidate  unique  millimeter  wave  properties  of  the  atmosphere. 
However,  the  model  is  by  no  means  complete  and  work  should  continue  to  perfect  it 
towards  more  specific  telecommunications  requirements.  We  recognize  a  need  for 
adding  subroutines  for  the  following  effects:  (a)  wet  snow  attenuation,  (b)  dust 
and  sandstorm  attenuation,  (c)  dispersive  rain  delay,  (d)  spectral  line  data  for 
trace  gases  and  man-made  pollutants,  and  (e)  statistical  data  on  spatial  and 
temporal  distributions  of  the  input  par  meters  over  scales  from  centimeters  and 
milliseconds  to  kilometers  and  years.  Implications  and  potential  of  the  MPM  for 
performance  evaluation  and  forecasting  require  more  attention,  and  integration  into 
system  models  is  suggested  to  allow  optimal  compromises  between  environmental  and 
technical  constraints.  These  efforts  should  proceed  hand-in-hand  with  extensive 
testing  and  updating  using  data  from  carefully  calibrated  field  tests  and  well 
controlled  laboratory  experiments. 
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APPENDIX  A. 

DETAILED  ATMOSPHERIC  (0  to  100  km)  ATTENUATION  AND  DISPERSION 
IN  FREQUENCY  RANGES  DOMINATED  BY  OXYGEN  LINES 
(52  to  67  GHz  and  117  to  121  GHz) 


The  catalog  of  graphical  results,  Figures  A1  to  A15,  is  /a  1  id  for  the  U.  S. 
Standard  Atmosphere  (1976)  assuming  a  relative  humidity  of  50  percent  between  0  and 
8  km  height.  Shown  are  responses  at  height  levels  h  =  0,  5,  10,  15,  20,  and  30  km 
obtained  with  program  PI  over  a  limited  frequency  range  (typically  400  MHz),  fol¬ 
lowed  by  Zeeman  patterns  of  0?  lines  which  lie  within  this  range.  The  three 
principal  attenuation  patterns  a-j(ir)  a2  3(0")  for  lines  identified  by  the  quantum 


1C  (Table  1)  are  calculated  over  the  height  range  h  =  30  to  100  km  using  program  P2. 
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Figure  A1  to  A15.  Attenuation  a  and  dispersion  D  over  15  frequency  ranges  between 
5?  and  67  GHz,  and  117  to  121  GHz  displaying  the  02  lines  1C  =  1  to  29  over 
the  height  range  h  =  0  to  30  km.  Also  shown  are  associated  Zeeman  attenuation 
patterns  for  altitudes  h  =  30  to  100  km. 

Note:  Each  frame  displays  tt(c^),  a+(a2),  and  a'(a^)  patterns  for  the  magnetic 
field  strengths  H  =  0.3  (left  hand)  and  H  =  0.6  G  (right  hand)  which  are 
symmetric  with  respect  to  the  center  axis  (interchange  o*  and  a").  The  fre¬ 
quency  deviation  Av  =  (f  -  v  )  is  between  ±40  MHz  for  H  =  30  and  ±2  MHz  for 
h  =  60  to  100  km.  The  maximum  attenuation  rate  aQ  is  that  of  the  isolated 

unsplit  line  (H  =  0);  the  value  in  parentheses  is  calculated  with  program  PI. 
The  values  of  a-j(v  )  in  decibels  per  kilometer  are  for  0.3  and  0.6  gauss. 
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Figure  A9.  continued 
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Figure  A10.  continued 
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APPENDIX  B. 


AB  INITIO  CALCULATION  OF  0£  LINE  STRENGTHS 

We  present  here  briefly  the  assumptions,  and  the  molecular  and  fundamental 
constants  that  are  used  in  an  ab  initio  calculation  of  line  strengths  (8)  for  O^ 
lines  in  air.  The  absorption  spectrum  N"(f)  of  a  pressure-broadened  line  is  given 
in  molecular  terms  by 

n  2i a~6u  o  9  -E./kT  -E/kT 

N"(f)  =  ^lFT  9£M2!4>|2(e  -e  )  F(f)  ppm  (Bl) 

For  spectral  lines  in  the  millimeter  wave  range,  the  approximation 

VkT  -VkT  -VkT 

e  -  e  =  ( hvQ/kT}e 

is  valid,  which  leads  via  (see  Equations  8,17) 


N"(f )  =  =  SF“ (f )  ppm 


to  the  definition  of  a  line  strength 

S  =  2irlO’6v0M  g£U2|$!2[exp(-E£/kT)]/3kTQ  Hz  (B2) 

16 

A  numerical  evaluation  of  (82)  for  the  linear  rotor  molecule  0^  (99. 76%)  in 
air  (20. 946%)  assumes  the  ideal  gas  law 

M  »  0.20946  *  0.9976  p/kT  =  1.5135  *  1019  p/T  cm'3 


for  the  number  density  of  molecular  absorbers;  and  the  appropriate  molecular 

parameters  are 

9*  58  1 

u2  ■  3.44034  x  io'40  (erg/Gauss)2 

6  •  43.100518  GHz 

E./kT  -  hBK(K  ♦  1 )/kT  «  2.06852  K(K  ♦  1 )/T 
Q  -  3kT/2hB  =  0.72516  T 

h  =  6.62618  x  io-27  erg*s 

k  =  1.38066  x  )0’16  erg/K 


With  these  values  Inserted  into  (B2),  one  obtains 


S  =  1.0892  x  10_4vo(p/T3)|$|2exp[-2.0685K(K  +  1)/T]  kHz  ,  (B3) 

when  the  line  center  frequency  vQ  is  in  GHz. 

The  matrix  elements  for  microwave  lines  are 

|4>+|2  =  kjK+(2K+  +  3)/K+  +  1) 

|*-!2  =  kj(K~  +  1 )(2K"  -  1  )/K~ 

and  k,  is  a  correction  factor  for  deviations  from  Hund's  case  "b"  as  listed  below 


K+ 

K- 

J 

kJ 

1- 

0 

1 . 0000 

1  + 

3- 

2 

.9808 

3+ 

5- 

4 

.9941 

5+ 

7- 

6 

.9972 

7+ 

9- 

8 

.9984 

9+ 

11- 

10 

.9989 

11+ 

13- 

12 

.9992 

>13+ 

>15- 

>14 

1.0000 

The  spectroscopic  coefficient  a^  is  defined  at  the  temperature  T  »  300  K  by  means 
of  (B3)  to  (Table  l) 

a^  =  S(300)/p  =  4.0342  x  10’1ZvQ|<»|2exp[-0.006895K(K  +  1)]  kHz/kPa,  (B4) 


and  Equation  (B3)  is  reduced  to  expression  (19),  which  is  used  for  the  Millimeter- 
Wave  Propagation  Model  (MPM).  These  involved  calculations  have  been  reported  many 
times  with  more  or  less  slightly  differing  constants  and  definitions.  This  appendix 
is  meant  to  reveal  the  extent  of  simplifications  and  reductions  that  had  to  be  per¬ 
formed  before  the  millimeter  wave  properties  of  dry  air  could  be  entered  into  an  MPM 
version.  The  expert  recognizes  our  numerical  assumptions  in  calculating  a1  and 
however,  we  refrain  from  defining  the  numerous  terms  (see  Section  2.1.1)  but  refer  to 
standard  references  (e.g.,  Waters,  1976;  Endo  and  Mizushima,  1982;  Poynter  and 
Pickett,  1981). 


APPENDIX  C. 

COMPARISON  OF  MODEL  (MPM)  PREDICTIONS  WITH 
SELECTED  LABORATORY/ FI ELD  DATA 

The  millimeter  wave  propagation  model  (MPM)  is  a  conglomerate  of  individual 
routines  addressing  formally  via  (7)  and  (8)  line  sums  of  02  and  H20,  continuum 
spectra  of  dry  air  and  H20,  as  well  as  fog  and  rain  contributions.  Each  routine 
has  been  developed  separately,  based  on  underlying  theory  and  supportive  experi¬ 
ments.  This  appendix  presents  several  examples  that  are  suitable  as  test  cases  for 
demonstrating  the  performance  of  the  MPM  in  actual  atmospheric  situations.  It  was 
difficult  to  find  corroborative  data  of  sufficient  quality  in  the  literature  that 
might  serve  such  a  role.  Reliability,  precision,  and  limited  scope  of  supporting 
meteorological  data  often  compromise  the  accuracy  of  results  deduced  from  field 
observations.  Generally,  controlled  laboratory  experiments  provide  more  serious 
tests  by  s  ..ulating  atmospheric  conditions  crucial  to  MPM  validations. 

We  selected  some  recent,  and  one  set  of  rather  classic,  attenuation  results 
from  laboratory  (Figure  Cl)  and  field  (Figures  C2  to  C5)  experiments  conducted  at 
frequencies  (21  ,  32,  96,  T « 0,  and  337  GHz)  that  are  particularly  sensitive  to  the 
new  water  vapur  continuum  (24).  Continuum  absorption  has  been  a  major  source  of 
uncertainty  in  past  modeling  efforts,  including  the  MPM.  The  comparison  between 
laboratory  data  and  MPM  predictions,  as  depicted  in  Figure  Cl,  indicates  good  agree¬ 
ment.  A  second  comparison  is  made  with  data  from  field  measurements  on  horizontal 
paths.  Figures  C2  and  C3  demonstrate  a  reasonable  fit  considering  the  scatter  in 
data  points  that  indicates  difficulties  plaguing  an  absolute  calibration.  For 
example,  the  72  data  points  in  Figure  C2  were  obtained  over  a  tnree-week  summer 
period  In  Boulder,  Colorado,  as  hourly  averages  taken  during  stable  weather  condi¬ 
tions  (afternoons).  The  vapor  pressure  e  was  deduced  from  microwave  refractivity 
readings  NQ(e,p,T)  at  the  receiver  site  by  applying  (4). 

The  last  example,  Figure  C4,  is  More  complex  since  it  involves  zenith  attenu¬ 
ation  Az  at  two  frequencies  (20.6  and  31.6  GHz),  measured  radiometrically  by  means 
of  atmospheric  brightness  Tg(ll)  from  three  locations  (hQ  =  0,  0.8,  and  1.6  km) 
in  the  United  States.  The  authors  successfully  fitted  straight  lines  (LF)  to  the 
water  vapor  dependences  as  summarised  in  Figure  C4.  The  experimental  conditions 
were  simulated  with  the  MPM  by  assuming  the  U.  S.  Standard  Atmosphere.  Adjusting 
the  dry  air  cut-offs  AQ  to  the  calculated  values (BF)  had  the  effect  of  slightly 
changing  the  water  vapor  slopes  (Av/V)z  within  data  scatter,  which  allowed  a 
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LABORATORY  DATA 
(Moist  Air,  1013  mb) 


t»pmti;es-ab(,T«3  rFTTsS  5fK"?^0hJraa5tf?Uatit?n  *i "° GH*  a?d  three 

/hi  ?qi  c  v  /nla  I,  ,  ic;  303.5  K  (Knight  and  Llewellyn-Jones,  1982), 
and  Auiier*  ,a9$e"yn-Jo"es  and  Knight,  1981)  and  (d)  at  31.8  GHz,'318  K  (Becker 
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ta _ m 

dB  dB/cm 

20.6  GHz 


LF 

0.049(10) 

0.199(7) 

BF 

0.175 

MF 

0.060 

0.167(2) 

rms  =  0.030dB 

31.6  GHz 

LF 

0.115(6) 

0.092(4) 

MF 

0.119 

0.089(3) 

rms  =  0.016dB 


20.6  GHz 


LF 

0.044(8) 

0.192(5) 

BF 

0.175 

MF 

0.048 

0.169(2) 

rms  =  0.01 7dB 

31.6  GHz 

LF 

0.103(4) 

0.076(2) 

MF 

0.100 

0.076(3) 

rms  =  0.01 7dB 

20.6  GHz 


LF 

0.076(7) 

0.173(6) 

BF 

0.045 

MF 

0.045 

0.170(2) 

rms  =  0.01 4dB 

31.6  GHz 

LF 

0.096(4) 

0.064(4) 

BF 

0.090 

MF 

0.085 

0.070(3) 

rms  =  0.008dB 


LF  =  linear  fit,  BF  =  biased  fit, 
MF  =  model  fit  (MPM) . 


Figure  C4.  MPM  predictions  for  radiometric  measurements  of  dry  air  zenith  attenua¬ 
tion  A0(dB)  and  water  vapor  zenith  attenuation  slope  Ay/V  (dB/cm)  from  ground  level 
observations  at  h0  =  0,  0.8,  and  1.6  km:  (a)  20.6  GHz  (D.  Hogg,  F.  Guiraud,  and 
E.  Westwater,  private  communication,  1983),  (b)  31.6  GHz  (Hogg  et  al . ,  1983). 
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satisfactory  fit  to  model  predictions.  Only  one  data  point  (A  at  20.6  GHz, 
hQ  =  1.6  km)  out  of  twelve  is  not  explained  by  the  simulation. 

A  more  complete  picture  of  modeling  the  data  given  in  Figure  C4  is  exhibited 
in  Figure  C5.  Dependences  on  total  pressure  P  (height  h)  and  four  values  for 
height-limited  (h  £  12  km)  water  vapor  contents  V  are  shown.  At  20.6  GHz,  about 
2/3  of  the  water  vapor  attenuation  Ay  comes  from  the  22.235  GHz  H^O  line;  at 
31.6  GHz,  the  picture  is  reversed  in  favor  of  the  water  vapor  continuum  (24).  In 
the  first  case,  line  pressure-broadening  'defeats  increases  in  the  slope  Ay/V  that 
are  caused  at  31.6  GHz  by  the  e^-term  of  (24).  Data  at  20.6  GHz  indicate  a 
possibility  that  the  22  GHz  line  strength  (b^  =  0.112,  Table  2)  is  a  little  bit  too 
low. 

In  summary,  good  agreement  between  MPM-predicted  attenuation  and  selected  data 
taken  under  a  variety  of  conditions  (that  is,  frequencies  between  21  and  337  GHz, 

v  =  0  to  50  g/m3,  V  =  0  to  3.5  cm,  T  =  274  to  318  K,  PQ  =  1013  to  820  mb)  lends 

credibility  to  the  millimeter  wave  propagation  model  MPM  presented  in  this  report. 
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